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SLOPE STABILITY ANALYSIS BY A NEW GRAPHICAL METHOD a _ 


SYNOPSIS. 

~The stability of slopes is commonly analyzed by the “method of slices — 

which involves certain simplifying assumptions and,as normally applied, leads" 

to conservative values of the factor of safety. A method of analysis developed 

by A. _W. . Bishop leads to more realistic values of the factor of safety but in- 


volves a tedious tabular form of computation. 


_ The derivation of Bishop’s equation for the factor of of safety of a slope is 3 out 


“a . This graphical method is 7 
then applied to the stability analysis of two hypothetical slopes.  ! eager 


| 
Notation.—The letter in this paper are defined where they 
first appear, in the illustrations | or in the > text, and are arranged alphabetical- 
for ‘convenience of reference, in the “Appendix. A 
General. —The determination of the stability of a heme is affected both Sl 
“inaccurate field data and by the inaccuracies inherent in the common methods 
of stability analysis. « Where the importance of the structure demands greater 
reliability from the field data, particularly pore- pressure m measurements, the 7 
use of greater refinements in calculation may be justified. 
- Note.—Discussion open until March 1, 1962. To extend the closing date one month, a | 
written request must be filed with the Executive Secretary, ASCE. This paper is part of — 
= copyrighted . Journal of the Soil Mechanics Division, Proceedings of the American — 
Society of Civil Engineers, , Vol. 87, No. SM 5, October, 1961. 
Senior in Civ. ‘Univ. of Adelaide, South Australia 
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It has been shown2 that in a _ typical earth dam, even assuming idealized 

elastic properties for the soil, local overstress will occur when the factor of 
safety (by the slip-circle method) lies below a value of about 1.8. _ As the ma- 7 . 
jority of stability problems involve factors of safety lower than this, a state of ¥ 
- plastic equilibrium must be considered to exist throughout some of the slope. cn " 

For these conditions the factor of safety (F) can be defined as the ratio of the 

available shear strength of the to that required to maintain 

when incipient failure is postulated.3 

If the strength mobilized is denoted as 5, 


which ct is the cohesion (in terms of the effective stress), gt denotes the ; an- ‘ 
gle of shearing resistance (in terms of effective stress), -p is the total normal 7 


stress, and uy repr esents the pore pressure. 


Sah! Failure is assumed along a continuous rupture surface, the shape of which 
-is influenced by the distribution of pore pressures and the variations ¢ of the = 


: shear parameters within the slope. For practical analysis a circular arcis _ 
_ usually nanan for the surface of of rupture and the problem is assumed to be 


_ DERIVATION OF OF BISHOP'S EC 3 EQUATION 
_ A complete graphical analysis of the equilibrium « of the. mass of soil above - 
‘the slip surface i is most laborious; several simplified procedures” have been | 7 
_ developed using the “friction circle” method which is, howev er, limited to— 
cases in which is constant over the whole failure surface. "More generally 
the “method of slices” is used with a simplifying simmation. about the effect | 
of the forces between the slices. 
The significance of this assumption may be examined by ‘considering the 
- equilibrium of the mass of soil (of unit thickness) bounded by the circular arc, 7 
_ ABCD, of radius R and with center atO [Fig.-1(a)}, 
aa the case inwhich no external forces act on the surface of the slope, a 
librium must exist between the weight of the soil above ABCDand the resultant | 


of the total forces acting on 


Enand the resultants of the total horizontal forces on the sections 


= the shear force acting on its base = sl; 
_ 2 The Stability of Earth Dams,” by A.W. Bishop, thesis presented to the Univ. of - 
London in London, England, in 1952, in partial fulfilment of the requirements: for the | 
of Doctor of Philosophy, po 
3 “The Use of the Slip Circle in the Stability oe ake of Slopes,” by A. W. Bishop, 7 


— 
— 
W the total weight of the slice of sol; 
the total normal force acting on its base =pl; 


SLOPE STABILITY 
= the breadth of the. element; 


= thea ngle between BC and | the 


x gre = the horizontal distance of the slice from the center of rate rota- 


2 becomes 
as 


(P - uy 1) tan 


Cc 


- onsidering the equilibrium of the soil in the slice above BC and resolving | 
ina direction normal to the surface [Fig. 1 (b) it it follows 


( Ep - 


Eq. 4 
Fe W cos a - uy 1)+ tan (Xn 


D. R. May’s of implies that the sum of the terms 
“representing the effect of the forces between the may be 
out serious loss of accuracy, that 


1+ ( 
sina“ 


— 
— q 
— 
— 
— 
| 
He 

= _ Equating th ht of the soil within ABCD with the _ 
moment of the Orme sliding surface we obtain 
— > 

| 
— 
— 
. — 
— 
— 
| May and 


FIG. FOROS THE SLICES METHOD 


Transformation 
for disfurbing 
forces 


RIVATION OF ‘STRESS AREA ORDIN. 


Transformed 
heighf of 


= 


Transformation 
for 


Bs 


PE STABILITY | 
In earth-dam design the construction pore pressures are often expressed 
asa function o of the wee of the column of soil above the e point considered, — 
@ J 


which B is a parameter on field data or laboratory tests 
and is assumed, for most design purposes, to have a constant value through - 
_ _ Thus, putting 1 = db sec @ the | expression for the factor of safety (Eq. 8) ) be- 
an W(cos - B sec «)] (10) 
Eqs. 8 and 10a are ‘implicit in accepted. methods of analyzing 
aoe stability. However, values of F obtained from these expressions are, in 
general, found to be conservative, especially where conditions permit deep ‘slip- 


circles around which the variation in @ is large. 
‘To derive a method of analysis that largely avoids this error, , the effective 
normal force (Pp - Uw 1) of Eqs. 4 is replaced by P' [Fig. 1 (b)] ‘and, after the 
forces on the slice have been resolved the expression is 


1- u cos a += sin a) 


Substituting Eq. 11 into Eq. 4 and 

= db sec 


db sec @ 


tan 9" tan @ 
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— before, , the term ini > 1) may be neglected for practical purposes 
without appreciable error and Eq. 13 then 
a | 
tan 
which is the form of Bishop’s equation 1 normally us used in the ar analysis of unsub- 

Eq. ‘11 indicates the fundamental difference between Bishop’s method and 
the original “method of slices.” Whereas the commonly accepted method ‘spe- 
cifically implies that the resultant effect of the vertical and horizontal — 
between slices is negligible, by resolving forces vertically | to obtain P' , Bishop _ 

_ eliminates the need for any assumption about the horizontal forces : and, for 
practical purposes, only the resultant effect of the vertical forces between 
slices is assumed to be negligible. 
_ The expression for F includes the value of the factor of safety on both sides 
of the equality, necessitating a solution bytrial and errorand requiringa time-. 

consuming tabular form of computation in the same as for the original 
method of slices. A purely graphical solution similar to that introduced by 
May for the method of slices ‘jiiiiaaes reduces alias tediousness | of the op- 


OF GRAPHICAL SOLUTION 


is conv to 
‘9 


+4 


: e in which form K can be | plotted for all generally used values of a, 9, ,» and F 


‘ale ‘unit of 
* the unit of rockfill: material; 


; oo = the height « of an elementary : slice; 

he = the height ¢ of the core material in an elementary slice; 


"f= the height of the filter material ‘in an elementary slice; q 
hy = the height of the rockfill material in an elementary slice; 
= - elementary slice for distrubing forces; 
= the | transformed height of the elementary slice for restoring forces; 


— 
4 
| 
— 
— 
— 
the unit weight of material, [Fig. 2 (a)]; 
= the unit weight of core material : - 
— 
db = the width of the elementary slice. — 
— 
= 


SLOPE STABILITY 
For ‘unsubmerged, homoge neous slopes, considering an 
[Fig. 2 (a)], expression for the factor Eq. 14, be > written 


ab h (1 - B) dbt 


: im zoned slopes [Fig. 2 (b)] it is convenient to handle the computation by 
_referring all unit weights to one basis, such as the unit weight of the — 
material (ve) The total the slice can, be written as 


in which hap is the transformed height of an element for disturbing forces | 


For Restoring pore in the restoring 
forces, the expression W (1 - B) becomes 


db he; 
hrp | 
y 
in which hyp R is | the transformed height of an for ‘restoring forces 
Thus, the process for various elements, full section can be 


into equivalent sections for the of the disturbing © 


and restoring forces. nee, 
sec a 


Eq. 16 then becomes 
\ 


db hyp sin 


The va values of c' and dg’: are those a “ini on the slip ci circle at ste partic- : 
assuming an initial value of F (conveniently F = 1.0), and, hence, estab- 
7 lishing an initial value of K, a simple graphical solution of the expression for 
the factor of can To do this, a cohesion circle is plotted 


— 
— 
he db + hg db hy db 
— 
RS 
— 


distance) vertically below the center of the slip (Figs. 3. 7, 
: 9). The cohesion distance will vary for the different materials through which s 
the slip circle passes and will, of course, be zero where the slip circle passes | 


through cohesionless material. ‘The cohesion | circle will then coincide with the © 


_ The basic graphical construction is shown in Fig. to which the following 


Let BA | denote the transformed height of an element for the case = 


Fras ale 
e 


els 


Surface 


gh 


“er 


He 


‘BE is plotted equal to AD. 


successive ve points on they whole section will give an area of stress 


between this line and the slip circle, representing, to scale, the > disturbing - 


forces on the slip circle (Figs. 6 and 8); that iy 
men ot of stress Sa db hyp sina .. 
~ Next, ain restoring ng forces, AF i is drawn at an in angle. or to BA to meet 
‘the horizontal through B at then rotated the BG. 


= hp tan ot (22) 


From Fig. 5 the value of corresponding to the values of and 


_ Pe 1. 0, is found and the height CG is reduced or increased in proportion (with 


— 
4.-PARTIALLY SUBMERGED SLOPES — 
= rom A a nernendicular is drawn to the slip circle radius BD and then 7 Pie | 
— 
i 
4 
— 


= 


: ith 


AN 


Wa 


am 


Yj 
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TALL AA iil | | 


5.- —CHART FOR THE DETERMINATION OF To the value of 
for @= 34°, ¢'= = 25°, and F = .. 5: _ Enter the chart on the left hand side at 
= 34°, proceed horizontally to intersect the curve for ¢'= 25°, rise 
to intersect the line for F =1.5 then continue horizontally to intersect the 
curve for positive values of ot and read off the value of K = 0.83 on the scale 
below. For values of @< 45° the factor K can be determined directly with 


equal facility on certain types of slide rule. 2. ee a 


SLOPE STABILITY ¥ &g 
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the use of proportional dividers or bya ‘simple construction 


is then rotated to the 


one = sec a K 


As before, joining successive points N on the whole section will give an 


_ areaof stress between the line and the cohesion circle representing th the restor-— 


_ sure coefficient for the distrubing forces (hp h =) and one including the coeffi- 


forces on circle ‘Tand 9), that 


of stress at db + db hpR tan tan K 


sit must be noted (Figs. 6 and 7) that, where the ile — passes through — 


- the core, there will be two transformed sections, one excluding the pore pres- — 


cient for the restoring forces (hpR).. Also, the restoring force stress diagram 


_ has discontinuities at points at which the slip circle passes from one material | 
_ to another (Fig. 7). This is because the value of tan 9° and, rr also, 
_«K K differs on either side of the vertical through the the point. 7 


_ (Fig. 3) in the ratio Kp /K, 9 in which Ky 9 is the computed value of K at the - 


The of safety is then determined 


Area of stress for restoring forces 


In fai this will ‘not, of course, be equal to the initially chosen value of 
=1.0 and new values of K must be determined, based on the value of F found 


to give a corrected area of restoring forces for use in a further determination 
of F. To do this it is necessary only to increase or decrease each vertical CN 
vertical with F = 1.0 and Kp is the revised value of K computed from the as- 
sumed value of F based on the first trial value. _ Two, or possibly three, trials” 
should suffice to obtain approximate coincidence between the chosen and com- — 


Considering the elementary slice shown in Fig. 4 (a) divided into two sec- 
tions by the free water surface e, a similar treatment leads to the expression © : 


for the factor of safetyas 


) + tan Wy +I 


7 
cOnesion Circle radius at M and 
— 
| im 
2 
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LOPE STABIL 


Transformed Seclion — 


‘FIG. 6.- -DISTURBING FORCES FOR A ZONED SLOPE (EXAMPLE 1) 


“Section 
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Te 
Centres of Cohesion Circles 


Ls 
of Slip Grele 


Scale (feet) 


SLOPE (EXAMPLE 2) 


Cohesion + of Cohesion Circle 


aler 


FIG, 9.- ~RESTORING FORCES FOR A PARTIALLY ‘SUBMERGED "HOMOGENEOUS 


SLOPE (EXAMPLE 
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— 
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STABILITY | 


“9 the of soil the slice the free water ace; 
= the pore pressure expressed as an excess over the 


= the unit weight of water; 
- the water elevation at the element; and 


depth of the slice below the free water er surface. 


it should be noted that the free water surface implies nothing about the mag- 
“nitude of the pore pressure inside the slope and is only used to obtain the sta- 
tically equivalent disturbing moment. The pore pressure at the base of the el- 
ement is still represented | by the standpipe water — it hy, as shown in _ 


Referring to Figs. 4 (a) and 4 (b), 


yy! = the _= the submerged unit weight of Serene, ‘below the free water surface; 


Yc = the sumberged unit weight of the core material; and 


= = the sumberged unit weight of the filter material. 


=y db (h 


=y, hy d 


ad) 
7) 
— 
i (he + hy - Zz) 
+76 Ne db+ 7}, (Z - he) db 
hy db + 744 (he + hg - Z) dl 


ot 
4 


in which hpp and hyp are the transformed heights of the element as before. 

From Fig. 4 and the expressions derived previously, the method of determin- 

_ TABLE 1,—PHYSICAL DATA FOR EXAMPLE 1 (FIGS.6 AND7) 


in 


900 


2,—RESU ULTSOF ANALY ES 

: 

Example 1 (Figs. Cand Factor of Safety 

Graphical Method 

Tabular Computation 


Example 2 (Figs. 8 and 9) 
Graphical 


ing the transformed heights hpp and is self- evident. clarity, o only 
hTR is shown in Fig. 


— 
— 
* 
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SLOPE STABILITY 


ms 
ad the d 
a. 28 thus re 


— db db hyp tan > K sec a 


* which is identical to Eq. 19 and can be solved in the same same way. 


ie slopes were selected and the factor of safety of pope: was calculated or 


Zoned fill no and B= =0.4 4 forth the core material (Figs. 6 and 

at 2) Homogeneous slope with 25 ft tailwater and with standpipe \ — eleva - 

tions represented by the broken line shown (Figs. 8 and 


Assumed 5 physical data for | example (1) are shown in Table 1 and for exam- _ 


ir ple (2) are shown in Fig. 8. Results of analyses are shown in Table 2. Table 2. 


CONCLUSIONS” 


ame graphical ‘method outlined is less tedious than the tabular computation 


The following symbols, adopted for use in the pa paper, ‘essentially 
with ‘ “Glossary of Terms and Definitions in Soil Mechanics,” prepared by the 
Committee on Glossary of Terms and Definitions in Soil Mechanics of the Soil 
Mechanics and Foundations depen Proceedings Paper 1826, , October 1958: 


B = pore pressure coefficient; eateiainenieain 


per unit area; 


= width of an elementary slice; 


= En = — on an elementary slice; 


a 
— 
— 
i 
. | 
— tm 
“Shes 
iy 7 - ual terms of the equation for the factor of safety. Once the transformed sec- ~ ae 
tion has been plotted the sequence of operations is almost as straghtforward 
1 as in May’s graphical method but only requires the measurement of one area © a 
— — 
 &§ 


— of an elementary slice: 


Pewee: of core material in an elementary slice; a 


= of filter in an | elementary slice; 


hy = height of rockfill ‘material in an elementary slice; 


height of standpipe water elevation | above e slip. circle; 


transformed height of an n elementary for forces; 

= = transformed of an elementary slice for restoring forces; 

K =; =a function of F, oe and 


L = height of filter’ material above free water surface; 


standpipe water. elevation n above free water surface; ae 
= total normal force e acting on - base of an eee alice; wet wn 
p : total normal stress; ie) Dek, 
R _ = radius of the Slip circle; 


= shear 


shear force acting on the base of an elementary 

= shear strength mobilized; 


ug * = pore pressure expressed as an excess over the hydrostatic pressure 
corresponding to the water level outside the slope; atti 
7s = = height of free water surface above base of filter material; 
We 


v total weight of an an elementary slice of soil: —— > 


= 


= weight of the portion of an elementary slice below the free water 

‘surface; 


= weight of ‘the portion of an in elementary, slice above the. free water 


= shear on an n ele 


uni walt weight material: 

unit weight of core material; 


- it weight o of material; 


= unit weight of ‘rockfill ‘material; 


— 
— 
— 
— 
— 
| 
| 
— 
| 
cetothe center 
— 
= 


- submerged unit weight of core material; z 


unit of filter material; 


— 
—_ 
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loading or column spacing. The values, although developed for continuous beam = 
footings, can also be applied to grid or mat foundations. The results obtained 


by: the empirical method are compared with those obtained by elastic analysis — 


and are found sufficiently accurate for practical purposes. | 
An expression for the rigidity of the foundation with consideration of | the 

“subgrade reaction of the soil is also developed. ‘This expression can be used, a 
within the stipulated limitations, as a characteristic for all foundations re-— 
-Sponsive to a contact pressure distribution governed by ‘subgrade | reactions, | 


design of combined footings of subgrade reactions has 


= stipulated regarding continuity, rigidity, ond variations in 


was the other reason stemmed from the difficulties of the theory and its 
application to practical cases, 


ae A considerable amount of research work was undertaken to explore | the pro- 


blems in connection 1 with the subgrade reaction produced by various soils. A 


_ Note.—Discussion open until March 1, 1962. To extend the closing date one month, a 
written request must be filed with the Ececutive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the American 
Society of Civil Engineers, Vol. 87, No.SM 5, October, 1961. ee | 
Assoc. ., Albert Kahn Associated Architects and Inc. ., Detroit, Mich. 
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? 
design pt purposes, ‘the of the subgrade reaction from the data 
_ tained from a standard soil exploration (see Appendix I). Because the com- 7 
_putations of moments,-contact pressures, and so on, by means of the exact 
- theory are not | sensitive to the selection of the k-value, ‘such an estimate is 
usually sufficient. If necessary, it can be checked later by actual testing. 
_ Much work ‘lie , done also to develop the basic idea theoretically. Since its 
original concept (Winkler, 1867), engineers were always intrigued bythe math- _ 
ematical potentialities of this method. Realizing its difficulties, tables 
been worked out (Hayashi, Hetenyi, and so forth) and methods have been de- : 
veloped (Malter, Gazis, Ray, Popov, and so on, to mention only a few) to fa- 
cilitate the "computations. . Most of the ‘papers produced quite ingenious new 


7 methods of approach to the exact sobation, but none appeared to be simple A 
enough to be used for general application. Consequently, most of the econom- 
_ ical advantage of the theory remained untapped. ae : 
_ Practical engineers realize that the uncertainty of their a assumptions often 
does not warrant the extreme exactness of the result. It is! known for example, 
7 soil in general does not satisfy the basic assumption of elastic support; 
that the k-factor is by no means a constant; that the homogeneity of the soil is 
certainly questionable, and that the rigidity of the super-structure is difficult 
‘, evaluate. These are only a few of the deficiencies that can be enumerated. 
It must be borne in mind, therefore, that the importance of this ‘method lies 
primarily in its value as a “tool, but not in the extreme exactness of its re- 


The following derivations will show that, if certain requirements are re ful 


filled, a simplified approximate method can be developed for the computation _ 
_ of moments, shears, and contact pressures for a combined footing which will 
- furnish . sufficiently accurate results. Other designs that cannot be fitted into 
the limitations and requirements of this simplified method must be solved by. 
- of the exact methods. This approach is similar to that of the American 
Concrete Institute (ACI) Building Code for flat slabs, which stipulates the re- - 
quirements of an empirical method for typical conditions alongside an exact 
_ __ Although the method described is for a continuous beam foundation, other — 
_ “combined footings, such as grid foundations and mat foundations, can be ana- — 
= ioeed on the same basis by dividing them into strips. Each strip, then, must — 
be considered as an independent unit and ‘must be analyzed by | using the full 
- column loads in each direction. In the case of a mat foundation, each strip will 
have a width extending from center to center of adjacent bays and will consist 
ofa column- “strip bordered by two half middle- -strips. The moments, shears, 7 
and contact pressures obtained by this method are for the entire strip. “They 
- will have tobe subdivided and apportioned tothe column- ~strip and toth the middle- 
is 


0 _ Basic Requirements. —The basic requirements for the applicability c of thi 


_gimplified method are asfollows: = 
>: 1, Each footing strip shall have not less than three bays or four ote - 


in a row in n order to ascertain the characteristics of a continuous beam. a 


Zz 3 “Evaluation of Coefficients « of Subgrade Reactions,” which was published by the a 


Institution of Civil Engineers, London S.W.I. 
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= COMBINED FOOTINGS 

a 2. The ratio between a adjacent spans ‘and adjacent column loads m may ay differ. re 
_ by a maximum of 20% of the greater value. It will be shownthat even such lib- 
_ eral deviations from the basic case of equal loads and equal spacings canbe 


tolerated without excessive loss in the accuracy of the end results. ae a 


a. The average column spacing 1 shall be be greater than | 1.75 — and nile 


tha: an —— in which is the “ characteristic” in the of a 


supported by subgrade reactions. T The meaning of this ‘requirement will be be ex- 
_ Ina continuous beam footing, ‘the column spacings must stay between cer- 
tain limits tokeepthe footing responsive tothe reactions of the subgrade . Col-— 
umn spacings that are too small will make the foundation act as infinitely | stiff, 
whereas column spacings that are too big will dissolve the continuous footing 7 
into single units with wasted interconnecting pieces, The selection of — 


4 


limitations permitted the computation of the maximum negative moment under 
the columns by means of an approximate formula using a simple straight-line 
_ expression instead of the trigonometric function of the exact theory, a 
= evaluation of the contact pressures, positive moments, and shears ay 
- simplified by the - assumption o of a a straight- -line distribution of the contact pres- 
- gure witha maximum value under the columns anda minimum value at the cen- 
terofeachbay, 
ee _ All formulae and derivations are based on the assumption that the continu uous 
beam footings o or strips are of prismatic cross section. Small variations in 
the width of the footing, making the narrowest part not less than three-fourths 
the widest part, will be tolerated without difficulty. 
Fig. 1 shows, schematically, the moments and the simplified distribution of 
the contact pressures on a continuous beam footing. But it should be kept = 


: mind th that the moments shown in ~~ are —— by the column loads only 


- 


— 
— 
| 
ee — 
(NY Mi | Mme 
= 
+ 
— 
a 
| 


‘and must be superposed, ‘if necessary, y, with moments caused by the settlement _ 
deformation (dishing) of the combined footing asa whole. 
= ‘This paper deals only with the moments and contact pressures produced by 


the 


__In order not to expand the scope of this paper unnecessarily, it will be pre- 
sumed that the classical theory of abeam supported by subgrade reactions (or, - 
as previously | called “beam or elastic foundation”) is known. Standard articles 
text books on this subject are available.4,5,6 
: _ Notation. —The letter symbols adopted for use in this paper are defined —— 
a, they first appear, , inthe illustrations or inthe text, and are listed alphabetically, 


for convenience of reference, in 


7 For a an 1 infinite long beam of prismatic cross-section, supported by ieeaiae 

reaction, the exact the following formula for the bending 


in which k isthe coefficient of subgrade ee of the soil, E denotes tl _ —" 
ulus of of the — to moment of in- 


4 “Theoriedes auf Elastischer Unterlage,* by K. Hayashi, Springer, Berlin, 


_ 5 “Beams on Elastic Foundation,” by M. . Hetenyi, Univ. of of ‘Michigan ‘Press, Ann 


" Die Lehre vonElastizitaet und Festikeit,” by E. Winkler, Prag., 1867. 


we 
{ 
— 
&g 
7 
qj inthis formula the “characteristic> 
— 


- and r for 1 various soil conditions and footing size ratios, see Appendix oe 


of AX. several loads are acting | on the beam, the total ‘moment a 
; _ point c will be the sum of all moment influences produced by the loads a ‘ 


* 
4X r 


. The total moment under an interior column load 


> 


“and the influence ¢ of one beprtens column load at each side, Py. 1 and Piet, lo- - 


om point a8 shown in in Fig 
spacing of the column loads. 
Cod concentrated loads are placed very close together, the foundation will ac act 
‘like an infinitely stiff beam, and the distribution of the contact pressure will 
: _ follow a straight line. In this case, the maximum negative moment at the col- 
will approach the value of the negative moment of an ‘inverted continuous 
beam carrying an equal distributed load (contact 


From Eq. ¢ the lower limit lmin can be evaluated as ——. For reasons of 

simplicity it was assumed to be ———. For column spacings closer than this 

value, ~~ 5 will furnish -moments that are greater than those produced by a 

_ straight- -line distribution of the contact pressure. Suchmoments are statically 

: _ produced only by loadings or contact pressures which are greater at the center 

_ of the span than at the support. However, , such a condition is not feasible in 
actual _— In addition to the selection of the minimum column spacing _ 


: with the maximum negative moment shall be checked so as not to exceed — 


4 this value is anes the foundation must be assumed infinitely stiff 
ane the of | the pressure will follow a line 


— 
> 
* 
=. 
X1, X2,... Xp are distances of loads Py, Pg... Py, from point c. iia 
If the spacing of the lo ds is hitrarily restricted for now to values he- 
6h(é<;7C73D tween and ——,, the t 
4 _ by a straight-line relation — 
| 
— 
| 
| 
— 
| 
il 
{ 
— 


a From the evaluation of the Qunetien for Cy x (see also Fig. 2), it canbe seen . 
that the moment influence from column loads located at a distance greater than 


= = can be disregarded. The was selected with —— to satisfy this 


requirement. At the same time, its magnitude eliminates the need for adding 


the: influence of the second _ adjacent column in either direction, even for i 
smallest admissible column spacing (3.5 = 2 x 1.75). 

Eq. 5 is an approximate formula, however, within the stipulated inccihahione 
it furnishes | quite accurate results, even if the permissible variations of col- 
umn loads and column spacings up to 20% of the greater value are considered. 


@ CONDITION 1 CONDITION 2 


_ loadings and | Spacings (as shown i in in Fig. 4) ys 


rT + 0.16) = - 
condition 1 ( Fig. 


A 
For 3 (Fig. 4c) 
| 


=- (1.25 Cal + 


— 
following equations give the relationship Detween he approximate q- 
— > 
@ CONDITIONS 
7 Ba 
(Ta) 
F3 ... (Tc) 


SM COMBINED FOOTINGS 
condition 4 (Fig. 


1 l 


TABLE 1, _—TABLE OF MOMENT FACTORS Fo TO F 


equal» column loads and | equal column spacings, the greatest deviation of the 
values found by means of the approximate Eq. 5 is about 3.5% of the exact value. 
_ The most unfavorable combination of loading, spacing, and stiffness ratio shows 


| a a deviation of 16%, but the average deviations are considerably smaller. At 


is the distribution of positive and — moments, but not the total amount, 
oy 


CONTACT PRESSURE UNDER INTERIOR COLUMNS 


_ The magnitude and distribution of the contact pressure under an interior 


column load is not determined directly from the exact theory, but fromthe mo- — 


ment obtained by Eq. 


In > the procedure, the following two assumptions 


1. —The resultant of all reactions acting on each in- 


= = (0.8 +Cygt 1. — 
1 6) \ ‘20.891 AO 
— 4 
i 
— 
1.75 | 0.5800 | 0.5960 | 0.6768 | 0.4950 | 0.6062 | 0.5914 | 0.6013 
| 0.6400 | 0.6414 | 0.7131 | 0.5517 | 0.6478 | 0.6282 | 0.6497 
sR 0.9400 | 0.9314 | 0.9451 | 0.9142 | 0.9205 | 0.9060 | 0.9310 E.* 
18 
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COMBINED 
(An interior foundation part is that portion of a combined which pri-| 
"marily supports one interior column load. It is assumed to extend from center 


_ Assumption No. 1 satisfies the fundamental requirement that the resultant 
of all subgrade reactions acting on a combined footing shall be equal, concen- 7 
gs and of opposite direction to the resultant of all column loads acting on the | _ 


combined footing; but, due to form, it simplifies the analysis. and 


facilitates the computation, — 


Assumption 2.—The pressures under an interior foundation part are 


_ assumed to have a maximum value under the column, a minimum value at the 
= of the bay, and shall follow a —— -line distribution between the col- — 
umn and center of the bay. (See Fig. 6.) — 
_ Because the negative moment at a column is affected but little by variations — 
- inthe loads of adjacent columns, the maximum contact pressure under acolumn > 
- can be approximated from a symmetrical loading condition over the average 
span l as _ shown in Fig. 7. From this: distribution of the contact pressure the 
negative moment | -Mj can be determined as 0 se 
Pm = | (3p Pi- + 5p 


; _ ‘Then minimum contact pressures at the middle of mor adja 


Pn ml * 1 Pi) 1) le (2 F + + 


= 


? 
4 | 
— 
ie 
— 
im 
— 
8 pj, the approximate maximum contact Pr 
| introducing Pm i he beam is obtained, as | va 
= ent spans can ther 
2p 
ok 


For equal s} spans and loads 

(18¢ 


By anaiihaaiiins each part of a combined footing, extending from center to 
center of bays, separately, two | minimum contact pressures are ob- 
tained at the middle of each bay. These contact pressures Pmr(1) belonging © ) 
to the left part and pmi(r) belonging to the right part are not necessarily alike 
because each was computed on a different basis. For the computation of the 
positive moment, they can be replaced by Bm as chown in Fig. 8. ._ 


~The moment Mo of the simple aquperine beam then becomes 


“and the positive moment 
16) 


in which Mj is | the ound verage of the | negative snenents under the columns é each 


| 
‘ar — final a the ‘the contact pressures, the shear can be de- 


_ Fig. 9 shows the. condition at the endof a continuous footing strip. The mo- 
s ment under the exterior column Pe consists of the moment due to the exterior | 
a column load itself, the influence of the next interior columns Pe+ 1 and the in- 
fluence of the free end. The first two parts can be treated | Similarly t to those 


for: an interior column. . The exact expression for the moment under an exter- 


ior columnis | = 


4 
— 
| 
— 
— 
= ag (Pili) +4 Pm Pi(r)) (18) 


‘FIG. 
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pression can be substituted and the e becomes 


The following equations give the relationship between the exact results for 


the various conditions of loadings, column spacing and st stiffness ‘ratios 1 in com- 
parison with those of the approximate Eq. 19, Pw 
Me (0.13 + 0.58) 
4X — 


For c condition 8 10a) 


For condition 9 (Fig. 10b). 


_ The e evaluation of loading conditions 8 to 10 for stiffness ratios \1 varying — oF 
esl 1. 75 to 3.50 is recorded in Table 2 and diagrammatically shown in Fig. 
__ 


T FACTORS F; TO F 


‘7875 ~+0,7980 0.7475 


0.8207 


0.9175 =| 0.9165 | 0.9332 | 
825 09825 | 09657 || 09725 
350 | 1.0150 | 0,9823 | 0,9858 0.9778 
The approximate formula furnishes results which deviate, for the most 
- wafavorable combination of loading, column Spacing, a and i stiffness, bya about 6. 6.5% 


from the exact values. 


90 SMS 

— 
(200) 
— 
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COMBINED FOOTINGS 


|, 


025 050 


— — 
4 
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int of a free end of a semi- -infinite beam supported 1 by sub-_ 
- grade reaction on the moment under the first column can be found through the — 7 
% application of a moment and a shear force at the free end, as indicated in Fig. 
_ 12, These two forces, called by Hetenyi “end-conditioning-forces,” make _ 
4 actual moment and shear of the infinite = beam vanish to simulate the con- 
dition of a free end. 


he influence of the end- nine forces « 2 the moment Me under the 


Cra 


. ‘Evaluating t this equation — the e help of tables given by Hetenyi the values ; 
ven in Table 3 and the | curve shown in Fig. 13 are obtained. — 7 


“Moment crosses the zero ) line at a distance 2 a= 8 n. If 


(0.1 13 Aly +1. 06 da - 


32. 
“hg | 
— — 
— 
— 
— 
— 
tions, moments which re- 
tor small end-projections, m ing. Subsequent 
especially h footing. Su q 
23 furnishes, rojecting part of the footing. Subse 
4 nigh subgrade reactions on the proj 


2 of the subgrade reactions for the exterior part towards the 
_ beam. _ This self- -balancing action led to the assumption that for small projec- 
_ tions the maximum contact pressure at the end of the projection (Pa) should not 
“exceed the contact pressure under the exterior column. This is shown in Fig. a 
om which indicates also the usually slight offset of the resultant of the contact 
pressures for the exterior footing part from the location of the exterior col- 
. To account for these adjustments, assumptions 1 and 2 have t been changed 


in 
to. assumptions la and 2a, as given below, if in with ‘moments 


contact pressures under exterior columns. 


Assumption la. —Under exterior columns, the of subgrade re- 
- actions acting on the exterior part of a combined beam footing must be equal nd 
opposite to the load of the exterior column. = 
: (An exterior foundation part is that portion of a combined footing which pri- 
marily supports an exterior column load. It is ascumed to extend from the 
center of the first bay y to the tip of the projection.) _ a eee 


! 
_EXTER. ‘FOUND. 


oe Assumption 2a. —The contact pressures | under an exterior foundation part 
-_" assumed to have a maximum value under the column, a minimum 1 value at 
_ the tip of the projection and at the center of the first bay, and shall follow a : 
-straight-1 line distribution in between. The contact pressure at the tip of the 


projection shall never exceed the contact pressure under the exterior column, 
adi The limitation for the contact pressure at the tip of the end projection leads: 


to a second formula for the negative moment under an exter ior column. From — 


Peat 


p 


2 Me = - Pe 2 4a+l 


— 

ee, 

max: Tp, Pe | aS DETERMINED FROM Pi 

q 
—. 

— 

——_ 


1961 
To assure the limitation the choice of Eq. 2 a 25 at shall a 


ased on la and 2a, the contact pressures ul under an exterior 


23 


wvern 


MUTUAL INFLUENCE OF ¢ ‘CONDITIONS 


aie be minimum length of a a beam, between exterior columns, is — 
can be considered as “long” ; 


beams and the end- conditions do not each other. 


s 
a 


* 
— 
+ Pa *Pm | ay > 
— 
— 
— ays, the 
selected for av s of clarity, 
on of the method. For s been omitted. 


‘COMBINED FOOTINGS 
‘Assume average 18 in. slab thickness (t), 0.225 kip pe per sq ft. Assume me soil 
to be stiff clay, allowable soil soil pressure is. 3000 ) psf., gi al aah tons per cu ft. ,. 


= 640 kips, Pe = 550 kips, 1 = 18 ft, b= 18 ft,a=5ft,k=2pkgi 


0 (see 1) ‘and fe= = 3000 psi. 


482 x 103 kips per sq ft, I= 18 x 1. 53/12 = = 6. 06 


—————— = 0.139 (see Appendix 1) a 
x 30 x 189 


‘1x A= 18 x 0.139 = 


_ 640 x 


2 Mi Sx 640 x 870 


= 2. 95 } —_ per sq a total unit it contact pressure, 


From Eq. 13¢ 


+0. 225: = 1. 45 kips per er Sq ft total unit contact pressure, - 


48 (Pitert * 


— 
lm 

— 

0.139 x 18 + 0.16) = - 870 ft kips~ 

7 = - Gy (0.2401 —— 

— 182 ( 49 + 4 x 22 + 49) = 1250 ft kips. 


= 
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Plan and Section 
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—~EXAMPLE ON MAT FOUNDATION 


ia 
— | 
tm 


al the middle bay an and 


COMBINED FOOTINGS 7 
Mm = Mo + Mj = 1250 - 870 = + 380 ft kips. 


@ 


(50 + 4x22 + 49) 


Design moments ot face of column 


FIG. 16.-SOIL PRESSURE AND BENDING MOMENT DIAGRAMS FOR MAT — ; 


or from Eq. 25 


CC 
— 5 
545 ft im 
| 
= 
#- Ta (0.13 1.06 da - 0.50) . - 


= 3.00 kip per sq per sq ft ft Lesaah unit contact pressure. 


| 


.’ a 4.2 kip per f 


| 
55 per sq ft total unit contact 


distribution a of the pressures and bending moments s is shown in 


_meenetes on Strap Foundation.—The following data is given (see Figs. 17 


and 18): Pj = 640 kips, Pe = 550 kips, 1 = 18 ft, a=5 ft, b= 6 ft 6in., t = 24 1 in. 
p= 0.67, ailowalile soil pressure = 8000 psf, Kgj = 7 250 (see Appendix I). 


4.147.200 


13.45 ft 25. 26.9 ft woe 


18 + 0. 16) =- = ft kips, 


= 


ft kipe for straight he contact pressure (886 < 960) 


= = 47 kips ft : 


+0. 3 = 55 per sq ft | total « contact pressure, 


a 
' 
0.13 
Mic 330 


A 


0" 
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= = 4.00 per sq ft total unit 


for the middle bay and 


4x 550 - 24x a). 


— 


1A, 


= 1280 886 = = 394 ft 


= 513 ft on 


x5+18 


= 48x 886 _ 3 x 640 _ 94 per ft . 
| 


6x 513 


Pe = 
and 

+0. 4s 8. 00 Kips F per sq ft total ‘unit ‘contact | pressure. 
690K 


Dead load excluded 


m moments a at face of of cokumn - 


Fise- Fa 
* 5.9 kips per sq ft total tal unit ¢ contact pressure, 
sown in 


? of 
_ 36.5 
— Pa = + 0.3 = 
_ The ea of the soil pressures | and id bending moments is ‘shown in Fig. 


6.5. 
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1S paper has shown that for continuous beam footings 
i be found byem- 
— tain limitations, the moments due to subgrade reactions can be found by ae BO 


making simplifying for the of the contact pres-_ 


sures under the footing, the magnitude of the contact pressures can be found - 
directly from the moments. The empirical formulae furnish results that are, 
for practical purposes, sufficiently close to those of an elastic analysis. il 


ings are not included in the formulae. 
_ Liberal deviations from equal column loads and equal column spacings can 


following 1 formulae and tables are based on work of Karl 
Hon, M. ASCE and Ralph B. Peck, F. ASCE. 7 od — 


+ 


value of p 


Side-Ratio. Factor p for on on Clay. = 


Stiff 


a b, “Evaluation of Coefficients of Subgrade Reaction” by K. Terzaghi, published by The = 
Institution of Civil Engineers, London, 1956. 
“Soil Mechanics in Practice” by K. Terzaghi and R. B. Peck, Wiley, 


_ Gridand mat foundations canbe treatedas cases of intersecting beam foun- — 
A ie - dations and can be solved by the same method. Examples illustrate the prac- 4 | De 
application of the method. = 
infeet| 10 | 15 | 20 | large] — 
TABLE 4.— COEFFICIENT OF SUBGRADE REACTION k 
— 
20-60 | 60-300 | 300-1000] 50-100 | 100-200 | >200 


values for the coefficient of reaction in Table 


2 


aa The following symbols have been adopted for use in in the paper and are ian aq 


length of end projection of combined 


compressive at 28 days, in ‘pounds per 


b, in feet to the fourth power; | 


Seen coefficient of vertical —_— reaction for a square area 1 al 
wide (see in tons cubic foot; 


= distance between adjacent columns, in 
of two adjacent column distances, in feet; 


moment at column for entire width of footing, in 


= 2 moment at interior column, for entire width of footing, in 


> 


average of column moments at each end Of a bay, in | Kip feet; 


: Mo = simple span beam moment (for the entire w width of the combined foot- 
- ing) for a beam supported by two 0 adjacent columns acted upon by con- ‘is 


— tact pressure due to column loads only, in kip feet; ae _ 
= number of blows in Standard acacia soil test; 


een = moment at center of bay, for | entire. width of combined footing, in kip 


— 
— 
— 
— 
— 
— 
average width of combined footing (b 
ar ters of bays or between center of 
bat 
N 


of exterior column, ng of weight of superstructure (D. L. 


= load of interior column, consisting of weight of eupetticemene (D.L. | 


+L.L.) excluding weight of footing and surcharge, in kips “a 


a L.L.) excluding weight of footing and surcharge, in kips; a 


oe contact pressure at tip of end projection due to column loads only, for 


entire width of combined footing, in kips per foot; hut ete pin 
‘contact pressure at exterior ‘column due to column loads o1 only, for en- 


tire: width of combined footing, in in kips per foot; 


_ tire width of combined footing, in kips per foot; 

ie contact pressure at center of bay, due to colume loads iui for en- rst 


verage the | Pm at the center of same as deter-_ 
mined | from adjacent t bays, in in kips per foot foot; 


total unit contact pressure, column loads plus footing weight 


= thickness of combined footing, in poreon 

characteristic in classical theory of beam by subgrade. re- 
ae, 


-size-or side-ratio factor of subgrade reaction Appendix I). 
wt imtad ha 
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‘DISPLACEMENT OF LO LONG FOOTINGS BY DYNAMIC LOADS 


oe A mathematical expression for the displacement of a long footing subjected a 

to a peak triangular force pulse is developed from theoretical considerations, 

_ The solution is presented in a manner convenient for the computation of a 


- Background.- _With the advent of of nuclear weapons, great interest has ne: 
shown in the study of construction capable | of protecting personnel and equip- 
ment from the effects of nuclear detonations. For the civil engineer the primary — » 
requirement is to design structural elements and structures capable of with- 
standing the high, time- dependent loads associated with the shock waves emitted . 
= such large energy releases. The rate of change of load within these blast — 
_ loadings is so rapid that the response of a structure would be appreciably 
different from its response to static loads of a similar magnitude. The — 
of such protective construction has therefore emphasized the study structural = 
_ Structural dynamics investigations had been pioneeredby s studies of ‘designs 
earthquake- resistant structures, and this work served as the natural start- 


ing place for further study. Both analytical and experimental work have a 


‘ ie: Note.—Discussion open until March 1, , 1962. To extend the closing date one month, a 
_ written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the American _ 
Society of Civil Engineers, Vol. 87, No.SM 5, October, 1961.0 
a. 1 Capt., U. S. Army, Corps of Engrs., U. S. Army Engr. Wa aterways Experiment Sta., 
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to determine the of structural materials, structural elements, 
a and total structures to blast-type loadings. Itis only natural that these studies _ 
pie were started with better-defined building materials such as steel and rein- | 
concrete, and left the less understood, nonhomogeneous materials such 
as soil to later investigators. Current knowledge of the dynamic response of © 7 
foundations is, therefore, greatly behind that of the superstructure. 
Purpose. and Scope. —Only | one aspect of the many problems involving the 
se properties and response of soils, that of the load-displacement re- 
sponse of long footings is presented, Specifically, the objective of the work 
described herein was to derive a mathematical expression for predicting the z ; 
permament of long footings in of soil 


Soils Division, WES, is invextigating the dynamic properties and response of 
soils through several research | studies being conducted or supervised by the 

-_* meetings involving representatives of the WES Soils Division and a 

: | =, various possibilities for mode of failure of footings under dynamic — 

: - loading were studied; among these were total shear failures involving rotation 
and punching of the footing, Rotation and punching as used herein distinguish | 
between failures where the failure surface is to one side of the footing (one- __ 

sided), causing the footing to rotate as it displaces downward, and where the 
failure surface is symmetrical about the center line of the footing so that no 

7 rotation occurs as the footing displaces downward. The writer felt that a study 
of a failure surface giving a punching-type failure would be worthwhile in 

order that the displacements predicted during rotational and punching failures 
_ under similar soil and loading conditions could be compared. Also, a failure 

7 surface for the punching failure could be chosen to allow eens com- 

for both frictional and cohesive soils. 


There was: reason to _ believe that a 
"prove more critical than the one-sided eetational surface, at least theoretical- 


“half of the two- sided surface were exactly similar to the one-sided surface, 
= mass of soil contained within the surface be twice as 


9 e 2 “The Behavior of Soils Under Dynamic Loadings,” by R. V. Whitman, Final Report 
_ of Lab. Studies, Contract DA-49-129-eng-227, Massachusetts Inst. of Tech., August, 1954. a 
em 3 “Studies of Investigations into the Dynamic Bearing Properties of Cohesionless — 
Soils,” by T. D. Landale, thesis, presented to the Massachusetts Inst. of Tech., Septem- 
1 4 “Design of Structures to Resist the Effects of Atomic. Weapons; Strength of Mate- 
_ vials and Structural Elements,” U. S. Army, Corps of Engrs., EM 1110-345-414, —_ 
_44 5 “Report on Survey of Literature in Connection with the Dynan ic Bearing Capacity 
of Soils,” by N. Kachaturian, Contract DA-22-079-eng-240, Univ. of Illinois, October, 


t... “Analytical Study of Dynamic Bearing Capacity of Found: ations,” by G. E. Trianda- 
Contract DA-22-079- Univ. of Illinois, June, 1960. 


_ writer became interested in the dynamic behavior of footings while on 
| 
| 
i 
eo i? failure surfaces can be shown to be quite small. On the other hand, the mass of © e Ae 
- 
— 
> 
— 


_ ther and therefore be the more 


mero 


hee General Approach.— The predicted permanent displacement (A4MAx) of a 
«Tong footing subjected to a vertical dynamic load is obtained by deriving an 
equation for footing displacement (4) versus time (t) and then solving the equa- - 
tion at the time of maximum displacement (tmMAx). The maximum occurs when 4 
_ the e velocity of the footing first seagsaperrsen zero. After this time, the motion | — P 


THIN FAILURE SURFA 


along the failure’ ‘surface. F Footing displacement. ceases at at this time because 


motion. As mentioned subsequently, this assumption does not consider elastic 


equation relating footing displacement | to time is obtained ned directly 
from the differential equation describing the n movement. The dynamic resisting is j 
forces s acting on the ‘soil mass and footing are expressed in differential form 
as effective resisting forces acting vertically at the footing. These effective 
Feaiating forces are equated to the looting function to to obtain the differ- 


: the reversed shear forces available are greater than forces acting to cae 


if 
— 
— 
+7. 
A 
— 
— 


he 
equation. The total on the is considered to comprise: the 
dynamic load, the ultimate static shear resistance of the footing (p), the inertial 


force of the accelerated soil mass (IF), andthe displacement-dependent restor- 


= force created by the movement of the center of gravity of the soil mass 
wn. This procedure will yield a second-order, nonhomogeneous differential 
equation for which the particular solution can be obtained by conventional 


Assumptions. —Solution of problems relating to the static bearing capacity 
of footings is complicated by the nonhomogeneous, anelastic properties of the — 
soil. In attempting to develop equations for the dynamic response of eee ; 4 


the problem is further strain- rate effects on shear strength, 


- Footing. —The footing is assumed long r relative to its width (2 B), ‘ain th 


pod effects are neglected. Although square and ‘rectangular footings a are of 


of of the and a brief analysis of each follow: | 


inherent i in the two-dimensional system. The footing is assumed to 
weightless and to impart a uniform loadtothe soil surface. The assumption of a 
_ weightlessness is conservative since the force required to accelerate the foot- 
would decrease the load on the soil surface. If deemed necessary, the mass" 
F of the footing and other superstructure could be included in the system, how- | 
ever, this is not within the scope of this paper. 
pee 2. Failure: ‘Surface. —Bearing failure occurs when the shear stress at all 
- points along the failure surface has exceeded the ultimate shear strength of | 
the soil. The dynamic failure surface assumed is the static failure surface = 
first described by Prandtl and later modified by K. Terzaghi, Hon, M. _ ASCE _ 
(Fig. 2). It is composed of two | spiral ‘curves of varying radius (R) starting at | - 
the point of a wedge-shaped zone located beneath the footing and terminating 
,) in two passive wedges which intersect the soil surface at an angle of 45°- ¢/2, 
in which | equals the angle of internal friction for the soil.7 8 As stated pre- 


viously, the primary reason for selecting this failure surface is that it hypo- rr 
thesizes the downward movement of the footing without rotation, Toe ay 
$3. Center of Spiral.—The center of each spiral surface (0) is assumed to 
lie on line ad (Fig. 2) or an extension thereof. The selection of the — 
os location of the spiral center is a prerequisite to precise computation of the | 
: ultimate shear resistance and of the mass of soil contained within the failure 
- surface. | The published data commonly used in determining the static bearing 
capacity are not based on an exact determination of the critical spiral center 
. ad because simplifying assumptions are available which do not greatly affect the 
_ accuracy of the static computations. The same assumptions are not believed 
pt permissible in the dynamic solution because ¢ of the added significance of the 
moving mass contained within the failure surface. The added labor to obtain — 
_ the more accurate failure surface and mass is believed quite necessary. a 
_4, Stress-strain Relation.—The soil is assumed to be a rigid-plastic ~ S 
i material, that is, the ultimate shear strength is reached without deformation _— 7 
a 7 ae Mechanics, "by I K. Terzaghi, John Wiley and Sons, Inc., New York, - 
8 “goi) Mechanics, Foundations, and Earth Structures, by G. G. 
MeGraw-Hill Inc., , New Y York, 1956, 221 
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ns and then the soil yields w with | no further increase in shear strength (Fig. o. 


The: discrepancy between the ; actual and computed deformations because of this — 
— should be small as long as the deformation is relatively large 


failure, This assumption implies that | the soil is incompressible and that 
4 5, Ultimate Shear Strength.—The ultimate shear strength is assumed equal 
= to cohesion (c) plus the product of normal stress on the failure surface and the 
tangent of the angle of internal friction. The ultimate | shear ‘strength is also” 


and must be incorporated into the solution, For nein. the ultimate al 
resistance could be multiplied by a strain-rate coefficient and substituted into 
the solution in the same manner, if the strain-rate effect is constant throughout 
the loading cycle. If the effect is more complicated, the incorporation |may 
require substantial changes in the differential equation. = 
Soil Mass and Movement.— The soil is assumed tobe of uniform density 
and to deform as a rigid body. This assumption is necessary in order to com-. - 
pute the inertial force, and ‘Should be approximately correct provided the dis- - 
__ placement is not excessive, that is does not destroy the initial geometry. 
Loading Function, .—The dynamic load applied to the footing is an 
bs triangular force pulse (Fig. 4). Suchafunction is convenient for steonate 


‘ions and approximates the probable footing loads under some practical condi- | 


mt affect the over-all procedure. | Another function could o oe » substituted, if it is 

felt that such would better represent the actual load applied to the Avabecegl 
- 8. Effect of Displacement on Ultimate Bearing Capacity.— Any downward | 

movement of the footing will displace the effective center of gravity of the soil 
=? _ contained within the failure surface in such a manner as to increase the ulti- 

_ ‘mate bearing capacity. This is not accounted for as a change in bearing capac- 

ity, but as adisplacement-dependent force acting to decrease the displacement, 

y that is to restore the footing to its original position. Displacement will also 
slightly shorten the length of shearing surface from the original ginal assumed ; 
failure. surface, but this effect is neglected in the computations, 


all ‘simplifying 2 assumptions are together, they certainly 
:- restesetions on the footing-soil system Geacribed. Nevertheless, it is be- 


that the sy ‘stem will accurately predict small displacements of footings. 
Cc. the predictions of long footing displacements under vertical dynamic — 
loads should be approximately correct if (1) the displacements are large com-_ 
- pared to displacements associated with static failure, and (2) the strength 
‘properties assigned to the soil closely describe the actual dynamic properties 


Notation. .—The for use in paper are 
here they first appear and are arranged alphabetically, for convenience eof 


reference, in the A endix, 


— 

— 

1e rate oI strain 

TOF Whe OF paper, tobe Independent OF tainl 
a> f strain-type effects is quite limited. Certainly, — 

since the present knowledge of str yp 
‘the rate-of-strain effect will modify the footing displace 

= 

= 

at 

— 
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In the following paragraphs mathematical for the forces 

7 on the soil mass are derived as effective forces acting vertically at the footing. _ 
These forces are then equated to the loading function to determine the differ- a 
ential | equation for the footing motion. Since the assumed failure surface is bs 
‘symmetrical about the « center line of the footing, only one-half of the footing and 7 
_associated failure surface is necessary to compute allforces. A unit length of 

1 footing is also assumed for all computations. 9 
‘Ultimate S Static Shear Resistance.—The ultimate static shear resistance of _ 

_ the footing may be defined as the footing load that will bring the soil beneath | - 

_ the footing into plastic equilibrium, that is, the condition in which the shear stress f 

“Sy along all of the failure surfaces is at the ultimate (plastic) value. To determine © 

"the ultimate static shear resistance, a a condition of ultimate shear | stress is 

_ assumed along all failure surfaces and the footing load required to maintain 

is footing load is the ultimate static shear resist- 


-FORCES ACTING ON MASS abe 


_ Consider soil mass abe (Fig. 5), and equate the sum of the vertical com-_ 
a ponents - of all static forces acting on the mass to zero. These forces include a 
’ the ultimate shear resistance, the weight of soil mass abe, the cohesion along ‘ 
_ line ab (Cp), , and the passive pressure force acting on face ab (Pp). The a 
_ ‘sive pressure force is considered to intersect t line ab at an angle of 90°-¢, — 
and therefore includes the frictional shear resist2nce along ab due to the nor= 


Vertical forces = = 


ng-1, BI 


| 
— 
— 
— 
lm 
— 
9 “Theoretical Study of the Visplacement OF Long Footings Dy Vynamic Loads,” U.S. 
A — — Army Engr. ‘Waterways Experiment Sta., Miscellaneous Paper No. 3-418, March, 1961. [i 


October, 1961 
Eq. 1 1 in terms of known parameters rameters except for the 
earth-pressure force (P ). To evaluate the force is separated 
into two components and Pf specifically defined so that their points of 
application along line are Force Pp is the passive pressure force 
considering the soil to be weightless and to have cohesion and a surcharge 
rs pressure at the surface. Force Pp is the passive pressure force considering ~ 
os the soil weight and neglecting the cohesion and surcharge (q). The forces are : 
evaluated separately by equating to zero the moments from all forces acting - - ie 
on mass abdf (Fig. 6). The sum of Pp + Pp then equals 
Fy _ Considering mass abdf to be in plastic oquilforiam, the following forces © 


contribute to the moment about point 0(M°): 
al Py passive force df neglecting weight of 


Py = passive force along df when c and q are 
= cohesive shear resistance ice along 


= cohesive : shear resistance along ab; 


weight of mass bOd; 
= = weight of mass ab0. 


ai Note in Fig. 6 that the resultant of the normal and frictional shear forces 


(N) along bd is through | point 0 and therefore does not enter the moment com- 2 


- Determine Pp: By assuming the soil weight tobe zero, forces Wj, W2, Ws, 
Pa and P" become zero. . P! will act through the mid-point of line ab. 


=0=M, * Mc, * 
por Pa 
R, 


= - 


‘in whichk is the ratio of a0/B, dimensionless; R,is the initial 7 


in feet; is of — at q, in feet; 45 + — all 


ae 
4 
— 
q 
— 
r= 
= 
— 
| 
— 
— 
hesion and surcharge to be zero, forces 
Pq, Cs, Cp, and Pp become zerc it of 
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_ Examination of the equation ‘for Pp and Py, shows them to be dependent 
the value of k, which is the ratio of length ab to B (Fig. 6), and therefore lo- 
cates the center of the spiral which is the center of rotation. Obviously, the — q 
sated value of the center of rotation is that which yields the minimum pas- 
sive earth-pressure force (Pp = Py + PY and therefore the minimum value _ 
ultimate static shear resistance. Since there is no convenient way to obtain and D, 
_ substitute this correct value and eliminate k from the equations, it must be 
determined t by trial and error for each set of problem parameters. 
Bin. The equation for the static shear resistance is obtained by direct hn 
"stitution of the sum Py + Pp into Eq. 1. This equation is then modified into a 
more convenient form (Eq. 4) by separation of terms into those dependent on 4 


Ts and q. The terms N Ne, and Nq designate coefficients of static shear = 


take 


B+ KBsind) +e (Ry +k B) -35kKB) cos ¢siny 
— 
— 
— 
— ain? gl 
isthe internal angle between B,, and line 0d; and id 
a (R,-7kB)cos¢siny = 
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R, (3 tan ¢ sin y - cos y) 
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Ry cos” (9 


(5a) 


ny 


acceleration c of alte mass due to the downward 
"movement of the footing will c cause se inertial forces: pep will resist such move~ 


‘individual soil mass and thereby dependent on the displacements. As the — 
first- order approximation of the motion, masses abe, acd, and ab0 are assumed _ 
to translate, and mass b0d to rotate about point 0 (Fig. 7). The prime (') - 
-ignates point locations after displacement. Note that it is kinematically impos- 
sible for there to be an abrupt transition from translation to rotation along 
lines Ob and Od, making the assumptions not exact but approximate. 
ad After” _ establishing the type of movement of each soil mass, it is necessary 
to | express the relative displacements. This is accomplished by assuming ef 
‘soil to be incompressible, from which it follows that equal volumes of each 
mass will be displaced by: aunit of the footing. For example, 
volume | aca' c' ' must equal volume aea' e', giving a relation of the relative dis. - 
placement of “masses abe and acd. From similar volume considerations it is" 
geen that the two edges es Roandt R2 (Fig. 6) of: of mass a not rotate the come 


_ 


resistance =| 
&g 
— 


- considered to be the average of the angular rotations computed for posers edge 
= The effective total inertial force is obtained by combining the inertial is 
forces acting on each separate mass using energy considerations. To be equiva- 
lent, the total inertial force must do the same total work when displaced down- 
_ ward as the am oe the work done by the individual ee forces as they are 


7 


= 


FIG. + 6.— 


S ACTING ON abdf 


a 

footing movement. rresistance 
displaced by the same namic inertial shear 


+k 


- Restoring . Force. .—Displacement of the soil mass within the failure surface c+ 
- due to downward movement of the footing will increase the restoring moment 7 


of the esteemed provided the angle changes donot become excessive. Only ay 
masses b0d and adf are assumed to contribute to the increase in restoring — 
moment. Mass cdf is excluded from computations because of its consideration ; 
7 as a passive wedge and the shortening of length cd with displacement make its _ 
effect more problematical. ‘The increase in restoring moment from the two 
“masses is converted to an effective restoring force acting on the footing to 
7. retard motion. The expression is given as Eq. 9. The term Np designate 
due to gravity restoring forces. 
— 


| cos ¥ + sin +R,” (sin 3tan ¢ cos 6, 


| 
Loading Function. —The fourth force acting on the footing is the loading 
function. For this treatise, the applied load is assumed to be an initial-peak — - 
Q is er force pulse which decays to zero at time tg (Fig. 4). The peak load : — 
is expressed in pressure units. Si _ the function is discontinuous at time 
ta, two equations are necessa 


r 0 <t <t ae Loading function = 
‘enwnees Equations .—The differential equations are set up by —_ the 
7 “four vertical forces to zero, . There must be be separate equations for before and 
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Displacement Equations .—The > solution the differential will 
yield equations of footing displacement versus time. The forms of ~ me 
ie solutions of Eq. 12(b) and Eq. 12(d) are found t to = ah: 


C, cos (Kt) + Cy sin (K t) + 


2 
; and » Co, and C, are 
integration. ‘The coefficients Cy and Cg are evaluated by the initial conditions. 
and C, are evaluated by the conditions of displacement and 


velocity at tq as defined by Eq. Solution and substitution of the coeffici- 


aa Coefficients. —The coefficients 
are dependent only on values of and k. Using magnitudes of from 0° to 


— ading function is defined inthat manner, : 
— 
— 
4 4 
[1 - cos (K t)] [sin (Kt) - Kt]..(14a) 
[1 - cos (K ta)| sin(Kt)-1.. (14b) 
q __NUMERICAL SOLUTION OF EQUATIONS 
— 
— 
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FIG. 8.—NONDIMENSIONAL MAXIMUM DISPLACEMENT 
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puter. The values obtained are given in Table 1 for every fifth degree, A com- — 
‘plete solution for each degree can be found elsewhere. eee ficial 


‘Maximum Displacement. —As previously stated, the maximum displacement 


z from Eq. 14(a) and Eq. 14(b) is the predicted permanent footing displacement, — 4 


| 


since downward motion ceases at the time of maximum displacement and re- - ’ 
bound is not considered. The electronic computer was used to determine the | 
maximum for ranges of Q/p of 1.10 to 3.00 


Footing Displacement under a Specific Load.— Footing displacement under 
a load may be computed in manner. Given the 


unit soil weight, in pounds per cubic foot; a 
cohesive shear strength, in pounds per square foot; 
angle of internal shear resistance, in degrees; 7 
| = surcharge, in pounds per square foot; 
Q= _ peak dynamic load, in pounds per square foot; and 


tq = duration of load pulse, inseconds. 
Determine (1) the ‘ultimate static shear resistance ©); (2) the values of NR and 


Np/Ny and and (3) the of 4max: | 


(1) To ‘determine | the ultimate static shear resistance (p), assume a value 
of k (which determines the failure surface) and determine coefficients N.,,, a 
and Ng from Table 1. Compute p for this k by substitution into Eq. 4, epeat 
— this procedure _ with different values of k until the minimum value of p is 
found. This value will be the ultimate static meer resistance, and the assoct- 

= k will define the correct center of rotation. = 

Note: 

value of “Only: these ranges 2 are shown in Table 1. 
(2) To determine the values of Np and./Np/Ny, use ¢ and the critical value 
_of k found in step (1) and 1 determine NR and4/NR/N] from Table 1. ai wie 
Todetermine the value of compute K (K=4 INR/NyB), Q/p, tq K, 
and Nry/p. Determine Nry/P)4max . from Fig. 8, using Q/p and tg K to enter 

the plot. Divide(NRy/p) 4max by the constant (Nrv/p) to determine the value 
-Effect of Load Parameters on Footing Displacement. —The effect of 
dynamic load or load duration on displacement canbe determined by repeating» ar 
procedure given previously for several values of or tg K. Plots of 


: pen gl versus Q or tq may then be — to ‘show the relation of 


and of k for the region where the ultimate static shear resistance could 
abe minimum the ents were evaluated on an IRM 650 dicit: 
— 
é 
— 
— 
é 
| 
— | 
— | 


TABLE 2.—SOLUTION OF EQ. 6(a) FOR NONDIMENSIONAL 


Nr Y 


0,300 0. 00004955, 
0,0005482 
0. 004748 


0.0001858 
0 001665 


a 


10,000 


0,01480 


1,873 
.003637 


0.04268 


0.3843 


1,505 


0.006227 


3.000. 2.215 
10.000 850 


ag 
Ff 


Solution of Sample Problems.— 
_ Problem 1.—Determine the maximum footing displacement for the following : 
footing, soil, and loading | Given: 
= 10 a= 100 Ib per sq ft 


60, 000 Ib per mas ft 


= 
ta 


Table 
for k = 0, = 13,274 + 25, 134 + 1,272 = 39, 680 
k= -0.10, ‘p= 11,354 + 25,553 + 1,292 = 38,199, 
c= 20, 9,887 + 27, 175 + 1,367 = 38,429 


= 38,052 lb per sq fi ft at k 0.15 


2. ‘Using Table 1 
= 6.9462 and In, /N 6.0748 
R 


. To max! 5 


J 
6.0748 0.50 x 1.921 = 0.96 


= 2. 2 ft 
Problem 2. —Determine the effect of various magnitudes of peak dynamic ae 


“load on the maximum footing for the soil, and 


ry 


ont wes 


For k = 0, p = 0 + 5712 = ADs - 


f k= 
Therefore p = Ib per ft at 


ie 1. p=N, By +N c+N, 400 Ny 


— 
— 
fa 
— 
— 
— =4ft aw 
Y = 100 lb per cu ty = 0.82 sec a 
q 
x 4 eae 


Note: For 0, there is only one — of tl the critical 


1.972 and 5. 588° 


. To determine A max as a function of Q: 


5. 588 _ 


pe 3. —TABULATION USING FIG. 8 


m 
| 

140 

| 8,996 
(11,424 
12652 


+ | | 


. 9.— FOOTING DISPLACEMENT—SAMPLE 20 
om ‘Using Fig. 8, , construct a tabulation as given in Table 3. eediec of sample 


problem 2 2 are plotted in Fig. 


COMMENTS: raw 


- The primary objective of of developing a mathematical expression for the 
displacement of long footings based on a punching failure and applicable 


— 
— 
ty K = 0.32 x 2.794 = 0,894 
— 
— ae — 
woe 
— 
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2 complex physical phenomenon, The assumptions have been restrictive to the 7 7 
- extent that precise correlation between the predicted and measured values of | 
_3 footing displacement should not be expected. At this time even the magnitude of 4 
_ the deviations cannot be assessed because of the lack of experimental data. 
P Laboratory tests are needed to furnish such information _ and to supplement 


AREA=53.8 FT? 
PERIOD =2.25 SEC 


/ FELLENIUS 
AREA =| 61.3 FT? 


-punching-type failures. 
using similar assumptions and of the displacements 


predicted in sample problem 2 with the predicted displacements for the same 
problem using a rotational failure shows the displacements for the punching 
‘ihemo to be larger.10 For example, at a peak dynamic load of 9000 lb per sq _ 
ft the modified Prandtl ‘failure surface gives a west matinee 5.5 times that for 


— 
the predicted displacements cannot be compared with measure g 
— 
"by G.E. Trianda- 
___ 10 “Analytical Study of Dynamic Bearing Capa 
| 


semicircular failure surface and 2. 3 me that Fellenius failure sur- 
a og = The shapes of these surfaces are shown in Fig. 10. Since the general 
on which the mathematical model is based will result in emailer 


‘the greater displacements associated with the punching 
that it is theoretically the more representative of the prototype. However, _ ] 
_ this does not mean that the punching failure using a modified Prandtl failure © 

surface is the most critical failure condition. i 
_ To determine why the punching failure yields larger displacements, a closer 7 

look at the three failure surfaces and the mathematical systems which repre- — 

sent them seems desirable. In the “Introduction” it was stated that the punch- — 

ing failure should be expected to give larger displacements, | since for similar oO 

failure surfaces the mass of soil displaced in the rotational failure would be 

larger. It is important here to realize that mass alone is not the only factor a 
affecting the response . of the system. With only a slight oversimplification it 

can be shown that two parameters, natural period (T) and ultimate shear — 

resistance, determine the response of the system. Assuming the ultimate shear - 
resistance to be equal, the system with the shorter period will have the larger 
displacements under the same load. The period is dependent on the resistance 


with displacement as well as the mass, and radius of ' gyration in wi in which the 


mass rotates wa i WR - Although the r mass is important in de- oa 

termining the period, it does appear be the predominant parameter. 


The area bounded by the three failure surtaces (Prandtl, semicircular, and 
Fe sllenius) and the natural period of each are ‘shown in Fig. 10. Even though 


"Prandtl, its period is closer to that of the semicircular failure, For the punch- 

ing failure, the soil mass is distributed closer to the centers of rotation, giving 
a smaller moment of inertia. Since the moment of inertia is in the numerator © 
of the a for the period, the effect of the two-sided failure surface is 


"SUMMARY 


to dynamic loads, assuming ultimate failure and 
-punching-type failure surface. The expressions are presented in a manner 
_ convenient for computation of displacement under initial-peak triangular load — 

; - pulses. Since the work is based on rather broad assumptions and is not sub- 
_ stantiated by experimentation, the expressions should be regarded as a re- _ 
_ search vehicle rather than adesign method for dynamically loaded foundations. > 


Predictions “of footing are compared for assumed rotational- 


that failure mode. “For both types were such that 
computed displacements should be theoretically less than the true displace- 


~ 
— 
wee 
— 
— 
. 
- 
| | 


DYNAMIC LO. 
ments, and therefore re it is concluded that the punching a 
approximation of wet true failure 
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‘Institute of Technology; R. B. Peck, F, ASCE, H, O. Ireland, M. ASCE, and 
GE. Triandafilidis, M. ASCE, University of Illinois; and Messrs. W. J. ‘Turn- : 
bull, ASCE, W. G. Shockley, F. ASCE, and R. W. Cunny, | ASCE, Water- 
ways Experiment Station, for their comments and advice. 


points used for r identification; 


_— half width of of footing, in feet; 7 
= of 


cohesive shear along line ab, in 


foot; 


distance from center line to point és passive wedge, in in 


Pie = acceleration due to gravity, infeet per square second; _ 


cleans I F = total inertial force, pounds per foot; 


= coefficient related to natural period footing; 


ratio of a0/B, dimensionless; 

= moment of forces | ce ,C,,, etc., about point 0, in poun nd-— 
= 


feet p per 


Ne 

= center of ‘spiral shear surface and center of rotation; 

ultimate static shear resi stance, in inpounds per square 


~ 


— 
— 
~ 
4 APPENDIX - NOTATIONS & 
al “Conesive snear strength of soil, in pounds per square 
— 
a Ney Ng N 
— 
— 


sree = passive earth- pressure force on line df when c and q | 
are zero, in per foot; 


per foot; 


passive earth-pressure force online ab soil 


are zero, in per foot; 


a surcharge on soil surface, in pounds per st square foot; 


= peak dynamic load, in pounds: per square foot; > 


radius of ‘spiral at in 


period, | n seconds; 
4 ty = = duration of load pulse, | in 1 seconds; 
tmax ax = tim time a at 


Wo, ete. = weight of soil masses, per foot foot; 
= = angle ab0; 


_ y = unit weight of soil, pounds per cubic foot; 


= maximum of footing, in feet; 


a0 A= downward transient displacement of footing, in feet; : 


line 


— 
— a 
= internal angle between Ry andvertical line through U; = 
62 = internal angle between Rog and 
angle of internal friction for soul, and 
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TANK FOUNDATIONS IN EASTERN VENEZUELA 


By Emery D. Ca D. Carlson! ; and Stephen P. Fricano, A. M. ASCE. 


_ The foundation experience obtained from ‘over 500 oil storage tanks con- 
structed in Eastern Venezuela is presented in summary form, The 20 yr settle- 


¢ 
“ment record of a twenty-eight-tank terminal constructed on soft clay soils is” 
& ven. It is observed that large steel storage tanks can be successfully operated | 
ow ith shell settlements as large i as 3 ft and bottom settlements as much as 2 ft 
- Continuous settlement records of three 150, 000 bbl tanks over a 10 yr. 
period, including controlled test loadings are shown, The three tanks are - 
_ located adjacent to each other on deep soft clay and have experienced unusually | 
large settlements. The foundation of one tank failed resulting in a ape oneal 
a ‘shell settlement of 5.75 ft with a differential settlement of 2.8 ft along the 


Shell. Soil conditions and analysis calculations of the failure ar are presented 


_ ‘Basic criteria used in the design of the company’s standard tank base is 
outlined. ‘Except in in the most unusual c nee a en earth base is used | 


Note. open until March 1, 1962. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part — 
of the copyrighted Journal of the Soil Mechanics and Foundations Division, Proceedings 
of the American Society of Civil Engineers, Vol. 87, No. SM 5, October, 1961. 
_ 1 Division Engr. of Eastern Venezuela, Creole Petroleum Corp., Caripito, Venezuela. 
Soils and Foundation Engr., Creole Petroleum Corp., Venezuela. 
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without “concrete ring wall topped by a 2 in. cold mix sand- 


_ From June 1928, when commercial oil deposits were discovered in Eastern — 
Venezuela, until mid-1960 the firm with which the writers are associated has | 
_ constructed tanks: ranging in size from 100 bbl to 150, 000 bbl capacity within 


4 San Joaquin 


Ciudad 


80 


= an area extending ateneeem from the Caribbean Coast to the Orinoco River 
>: and from the Eastern coast 180 miles inland. It includes the Venezuelan states 
of Monagas, Sucre and most of Anzoategui. The low mountains of the state of 
Sucre and the central plains area directly south make up a water shed area 
which drains toward the south and east to the Orinoco River, and toward the 
west to the Unare River. The eastern halfof this area, principally in the state 
nd Monagas, forms part of the great Orinoco alluvial delta which drains into 
the Gulf of Paria adjacent to Trinidad, andinto the Atlantic Ocean. The annual - 
; rainfall varies from 30 in. over the inland plain area to 80 in. to 100 in. in the 
Delta area, The The majority the tanks built built are located at Ciudad 


= 
— = SQ 
“ 
— 
— 
— 
— 


Barrancas, and -Caripito” (Fig. and are constructed on the 


Orinoco Delta sedimentary fine sands ,Silts, andclays.§ 
os 10st difficult soil conditions encountered in the Eastern Vene- _ 


zuelan tank farms are those inthe Caripitoterminal located next to the Caripe 7 
= two miles from where it joins the San Juan River. In this area ane 


aes data, obtained a as s early as 1940, has been ‘compiled for | approxi- 
mately twenty tanks inthe Caripitotankfarm,. 
The detailed time-settlement curves for the initial load tests and subse- 
any eight year service » records of three 150-ft diameter tanks are presented. ; 
The failure of one of the three foundations is explained with a description of — 
‘the remedial measures taken to restore it to service. The standard ‘compacted-— 
fill tank base, , used conditions permit, i is presented, 
SUMMARY OF TANKS CONSTRUCTED 
a re total of 555 tanks of all sizes are still in service in Eastern Venezuela, 
built by the company with which the writers are involved from 1928 to 1960, 
Numerous others have ‘been built and dismantled or replaced. Of the total 
present tankage, 395 are 1,000 bblto 150,000 bbl capacity, and 160 tanks range : 
in capacity from 100 bbl to 1,000 bbl. This tankage represents a total capacity 
of 11,173,200 bbl. Table 1 summarizes the tank sizes and types of foundations - 
built. It is particularly : significant that only twelve tanks were built on founda- 


=" 


on which they were b built. 


é 
tions other than compacted- fill oe one aie the varying soil conditions : 


20° YR SETTLEMENT ‘EXPERIENCE OF THE HE CARIPITO TANK FARM 


terminal located within a large curve of the Caripe River just two miles west 
of the San Juan River. This is shownin Figs . 2(a) and 2(b). Construction of the — 
terminal began in 1931. Additions have been made periodically until today 


there are 28 tanks with a total capacity of 2,752,000 bbl. The tanks range in | 


capacity from 30,000 bbl. to 150,000 bbl. with diameters from 70 ft to 150 ft. 
The dimensions , Capacities, , and elevations of each tank are shown in Table 2. 

Soil Profiles. —Numerous borings have been taken in the tank farm area 
since the 1930’ Ss. The soil profile shown in Fig. 3(a), is representative of the — 


underlying strata. The whole area is underlain with soft gray clay varying» 


from 6 ft to 30 ft in thickness. Under this there are intermixed sand and clay 
strata down to a depth of 90 ft where sand, ara, and boulders are en- 


“countered, 


Se ttlement Data.—It may be noted that _ the twenty tanks for which there | 


are elevation records in the 1940-1950 period, the average shell settlement _ 


0.46 ft. The same twenty tanks registered total settlements from 1940 to 1960 
from 0.00 ft to 1.24 ft with a mean of 0.66 ft as shown in Table 2, | ie 
All tanks were constructed on compacted earth bases except Tanks 22 and _ 
29, which have concrete ring walls. Between the years 1940 and 1960 these | 
two tanks had average shell settlements of 0.62 ft and 0.16 ft, ‘respectively. 
can be seen for the concrete ring that it t it pro- 


= from 0.03 ft to 1.16 ft. The mean settlement for the twenty tanks was pa 
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‘October, 1961 
- vides a level base for the original | erection of the tank shell. The r: ring may 
also provide some protection against differential settlements along the shell 


due to soft pockets. However, the experience in Venezuela has not revealed — 
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FIG. 2. _TERMINAL ‘TANK 

any etiaiaiaii tank shells resulting from —? ‘soft pockets along the shell of 
the other tanks constructed on compacted bases. The center portions of all 
tanks in this "the roof columns an 1.0 ft to” 2. 0 and nave 
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TANK FOUNDATIONS 


TABLE _SUMMARY OF TANKS BUILT IN EASTERN VENEZUELA 
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75x12 | 1,000 | 


Timber Piles 
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x45 ~ 139. 300, | 


Total Capacity—11,173,200 Bbl 
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October, 1961 
_ tinuous data on individual tanks is not available « except f for three 150,000 bbl 


SETTLEMENT RECORDS OF THREE 150,000 BBL TANKS 


The most interesting tank e experience in the Caripito Terminal has been 
the settlement performance of three 150,000 bbl tanks constructed in the south- 
east corner of the tank farm. Prior to construction of the tanks, soil samples | 
were taken fr from two deep borings at the site. Laboratory tests, performed in 7 
: the field, were limited by the test equipment available locally at that time . 
They inckuded blow counts of | moteture limits and 
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"classification. From the results o of these tests and from the general experience 
tank bases inthe same area, it was felt that settlement be 


with heights of 2.4 ft, 2.6 ft, and 3.7 ft for Tanks” 37, 38, and 39, respectively. 
dj The bases were constructed with a a sandy clay gravel mixture brought from 
a point 2 km (1.24 miles) inland. No record of compaction of base material is 


available. ’ The three tanks constructed were API "I (American Peti roleum | Institute) 


| — 
— 
ea] 
— 
+ 500 pst + — = 
sequently, both short and long term observation schedules were = 
record the movement of these 4 
Se Sle in in the general area slopes slightly to the southeast. In order _ 
— terrain in th elevations for each of the three tank 
to obtain approximately the same base e 


Average Shell 
Settlement, in feet 
Capacity, Diameter, 


in barrels | in feet 
| @® | @ 
55,000 | 116 | 04-32 
55,000 | 32 0.25 
__ 80, 000 42 35 
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October, 1961 
welded cone- roof type, 150 ft indiameter and 48 ft high. To facilitate settlement 

readings, clip angles were welded to the shells at the quarter points. a? a _ 

_ Tank No. 37.—This tank was test loaded with water beginning September 9, | 
- 1949, A total of 46.5 ft of water was pumped into the tank at a fairly uniform 
=o a 30 day period for a gross pressure on the sub-soil (including the 
tank at 75 psf, and the base at 250 psf) of 3,235 psf. The load was held for 90 
_ days from October 9, 1949 until January 7, , 1950. ‘The water was then n pumped — 
out at a nearly uniform rate in 16 days. 
The time-settlement curve in Fig. 4 shows fairly uniform settlement of the 
shell which averaged 60 cm (1.96 ft) after 90 days loading. Points 1, 2, 3 and 4 
_ are observation points on the tank shell. Their orientation is indicated. Maxi- 
mum differential settlement of the ‘shell was ( only 9 cm ©. 30 ft). Rebound of 
the shell after unloading averaged 5 cm (0.17 ft). Cross sections showing the 
settlement of the tank interior are shown in Fig. 5, with point orientation as 


shown in Fig. 4, Examination of the settlement — of the tank interior 
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shows settlements near _ center of the tank of over 50 cm m (1. 67 ft). more than 
“7 along the shell. No evidence of general shear failure is apparent, although some 
local shearing of the compacted base took place as the tank punched into it. 7 
Sg This tank was placed in service in 1950 and has, throughout the year, an 
: average | soil loading of 2,700 psf. The time- settlement curve from 1952 to 
< 1960 is shown in Fig. 6(a). Total additional settlement of the shell while in 
normal service averages approximately 42 cm (1.40 ft). 
aie Tank No. 38.—This tank was test loaded with water beginning November he 
1949. Again, approximately 46.0 ft of water was pumped into the tank at a uni- 
form rate ov er a 30 day period for a gross pressure of 3,230 psf. The load © 
was held for a a 60 day period from December 1, 1949 until February 1, 1950. - 
The water was then removed at a uniform rate in 15 days. The time- settle-_ 
- ment curve ‘Showed an an 1 increasing rate rate of se settlement as the point | of full full 1 load 
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was approached. iiiesaaion 20 days ys after the full load was was placed the rate. 
of settlement decreased and stayed so uniform that the test was discontinued 
after 60 days in order to place the tank in service. +The tank shell showed an 

average settlement of 82 cm @. 67 ft) at the end of the test with a maximum 
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the shell of 12 cm (0.40 ft). in the case of Tank 
_ No. . 37 the interior shows settlements of about 40 cm (1.33 ft) more than along 
the shell Fig. 6(a) . Again, no evidence of general shear failure is apparent © 


although the | tank into the base 30 cm cm (1.0 ft) 0 ft). In 1950 
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FOUNDATIONS 
it was s placed in | service similar to Tank 37. Its time- siiatiatiaas curve » from 
1952 to 1960 is shown in Fig. 6). Total additional si settlement of the ‘shell < 
Tank No. 39.—The tank was first loaded with water beginning Seeeetes 5, 
1950. A total of 45.75 ft of water was pumped into the tank at a uniform rate 
(with one 4 day interruption) over a period of 33 days for a gross pressure of 
3,350 psf. This load was held until February 2, 1951 when it was increased to 7 
3, 500 psf. and held for 82 days until April 25, 1951. The water was then re- 
leased by gravity flow. 
- ¥ _ An analysis of the loading and settlement curves I reveals a rapidly increas- — 
ing differential - settlement immediately after the 3,350 psf gross load was 2 
placed. However, after several days, the rate of settlement decreased con- 
. siderably until February 2 when the total load pressure was increased to 
— 3,500 psf. _ About four days following this load increase, the settlement rate of 
7 points 2 and 3 on the southeast side of the tank speeded up for a period of 7 
four weeks. The tank tilted badly to the southeast until it stabilized with a | 
mi maximum differential settlement along the shellof 85cm (2.8 ft). Total settle- 7 
ment on April 25 varied between 90 cm and 175 cm (3 ft and 5.75 ft) at the . shell, 
tank was considered unsafe for general service. = 
Examination of the settlement profiles. of the tank interior (Fig. 5) indicated 
greater’ settlements were concentrated along the tank edge between points 
2 and 3. Settlement along the southeast edge was nearly 1.0 cm (3.3 ft) more> 
_ than at the center of the tank. _A shear failure was definitely indicated. The — 
failure appeared» to be confined to the | base and ‘upper strata. However, no 
_ ANALYSIS OF TANK No 39 FOUNDATION — 
"Description n of Foundation Soils.—The original foundation of Tank } No, 1. 39. 
consisted of a relatively uncompacted fill of clayey gravel inthe shape ofa _ 
wedge varying in thickness from 2 ft to 4 ft. This compensated for gentle 
slope toward the southeast. The subsoil conditions are represented by the soil 
profile shown in Fig. 7. These data were collected from two deep borings taken © 
in January, 1951. ‘The | top 6 ft consists of a medium brown clay with undrained . 
‘shear values of about 720 psf. ‘Shear values are one-half the | compression — 
 vength from the unconfined compressiontest. Below this is a stiff blue ef 


yy down to - -19 ft with a an in undrained shear value of 1,860 pst. From - 19 ft ry 


‘From - -52 ft. down, shear strength3 increases up to 700 psf at the 
_ side of the tank. Penetration resistance required to drive a spoon sampler 12. 
in, with a 140 Ib hammer varied up to 14 blows. Water contents ranged oe 
30% to 74% of dry weight. Dry weight varied from 55 pef to90.8pef. 
_ Analysis of Failure Conditions.—The foundation analysis of the failure, : 


transverse shear tests.4 Overload ratios of the failure condition varied from 
“Notes on Basic Principles of Shear Strength,” by R. V. Whitman, 1959. 
4 “Design Memorandum on Bearing Capacity of Spread Footings on Conetive & Soil,” ; 
= 


4 45° ‘from ‘the horizontal, Shear values. were e those obtained from = 
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2. to 3. in the ‘These a: are approximately equivalent to 
Terzaghi-Peck®5 safety factors®4 of 1. 0 to2.0. The reconstructed tank founda- — 
tion was calculated to have an overload ratio of 2.5, , equivalent to a Terzaghi- 

i Because unconfined compression strength tests were made on the original — 
- samples at the same time as the Housel transverse shear values were ¥ 
taken, it is possible to make an analysis using undrained shear strengths 

equal to one-half the unconfined compression strength using the stability 4 
= formula of A.W. Skempton.? = 

_ The stability analysis was done by considering three cases. Case I assumes 
_a failure of the base and the soft top stratum from 0 ft to -6.0 ft, with the base 7 

having a shear value equal to the top stratum, 720 psf. Case II assumes a 

ss failure circle extending down through the next stratum to -19.0 ft, with a shear 7 

7 value of 1,860 psf. Case III assumed a failure circle extending down through — 

the soft stratum having a 


"Because the three assumed ‘failure cases varied from an annualar strip 
to. a full half -eireular: aise , the Ne value in the bearing capacity formula — 


as -proposed® by A. W. Skempton in 1951 was modified from 5.2 to5.69 as — 


suggested by L. Bjerrum. 8 The calculations are given in Appendix I for each 
case. The results are summarized in Table 3 
It is readily apparent from the safety factor ‘found in the failure condition 
tt that the foundation could have failed in both the 0 ft to -6.0 ft stratum and the 
_ =19.0 ft to 73.0 ft stratum. From field observations and study of the settlement 
_ Curves and sections of tank bottoms, as well as the stability analysis, the 


believe there three factors involved in the foundation failure. 


1. The tank the relatively poorly compacted clayey gravel 
"base, failing it and the top soft clay stratum . This was evident from field obser-_ 
a a 2. The foundation auiaun , -19. 0 ft to -73.0 ft, also failed in shear. oe! 
pe section showing the settlement of the tank interior plainly indicated the 
failure ap ee nearly to the tank center, approximately 70.0 ft from q 
al 3. The entire tank base settled considerably asaresult of rapid consolida- 
‘t.: of the underlying strata. The load-settlement curves indicate a fairly 
rapid settlement rate, howevera settlement rate versus load plot (Fig. 8) shows 
the settlement rates rapidly decreasing under constant load, It seems logical , 
to assume that the underlying strata were undergoing a constant consolidation 7] 
process with corresponding shear strength increases, oven immediately before 
and after the two periods of apparent overload and failure. 
Remedial -Measures.—The tank floated off the ‘base by flooding the 4 
area within the fire walls and then towed to a new location within the area. [ff 


= 
5 “Soil” Mechanics in ‘Engineering Practice, ” by K. Termes and R. B. Peck, John 
“Soil ‘Mechanics, Foundations, and Earth Tschebotarioff, 
McGraw-Hill, New York, 1951, (a) p. 405. 
“The ‘Bearing Capacity of Clays,” ” by A. WwW. Skempton, ‘Building Research Cong., 

= Foundation Failure of an Oil Tank in Feetetsheut, Norway,” by he Bjerrum and > 

Overland, 4th Internatl. Conf. on Soil ‘1957. 
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a iain the area the gravel fill was removed from the base and the cae 
brown clay layer from 0 ft to-6 ft was excavated. The base was then rebuilt to 
its original elevation with a well compacted shale. The area was then ao 7 
and the tank floated back to its original position. After the area was again | 
drained, an earth a varying from 4 ftto 6 ft was ‘placed around the tank © 
Although c: considerable long term settlements were expected in Tank No. 39, 
ig ‘due to deformation of the very soft deep clay stratum and consolidation of 7 
silty sandy clay stratum, it was believed that the rebuilt base and shaped sur- 
™ charge would reduce the differential settlement to serviceable proportions, — 
: The: safety factors, 2.98 to 3.02, obtained in the top strata by modifying the 
base, preclude failure in these strata. It seems that consolidation of the -19.0 ft 
to ) 8. 0 ft stratum, with subsequent increase in shear strength, ,isa logical | 
answer to why the tank foundation has performed satisfactorily after being 
modified . The safety factor of only 1.05 inthis stratum is based on the oxtafnel 
borings taken during the first load test and do not reflect the e subsequent cc con-— 


oe due to both the first and second load sees 


Checked, Safety | Bearing 
in psf | psf 


an 


a _ Load Test on the Modified Foundation. —The testing of the modified 
foundation began November 10, 1951 with the addition of 5 ft increments | of 


water. Each increment was added over a day pe period and left in place 
= before adding the next increment. This method of loading was used in 


2 


order to allow more time for consolidation to take place in the underlying | 
Strata than had been allowed in the previous test. It also made it simpler to — 
detect any changes in | rate of settlement which would indicate impending failure 
or unacceptable differential settlements. Actual loading rates and resulting © 
settlements are shown in Fig. 10. ~The total settlement of the shell ranged 
fom 15. 4 em to 21.4 cm (0.51. ft to 0.71 ft) with a maximum differential of — ; 
i: 0 cm to (0.20 ft). Since the settlements were so small, the full load of 3,500 zt 
psf was held for only 10 days. After the water load was removed, the tank re- © 
bounded an average of 3.0 ¢ cm (0.10 ft). Cross section of the settlement of the 

_ tank interior are shown in Fig. 11 . The long- term service time- settlement : 
curve, Fig. 6(c), shows an average additional shell settlement of 75 cm, 
(2.5 ft), and a maximum differential of | of 21 em 0. 7 ft). The 1e direction of tilting 
from ‘Southeast southwest. 
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TANK FOUNDATIONS | 
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factors contributing the successful performance of the ‘modified 
foundation are assumed to be as follows: 


Partial some sub-soil strata from the preloading action 


the two water tests; 
2. Replacement of the 6 ft. thick softbrownclay layer and the uncompacted 


clayey gravel fill with well compacted shale; 
3. Placement of the balanced surcharge around the tank. 
_ Although it was recognized that the addition of the surcharge rs in- 
it was believed to be justified in view of 


“crease total 


~Final elevation In after load test No. 2 


Fit nal elevation after oad te test No, 2 
SECTION 2-4 


Surface ‘primed andcoated 6in.to 12in.of top soil removed 
with RC-2 asphalt 


FIG. 12. .-—-TYPICAL ‘TANK FOUNDATION 

the fact that it would reduce the possibility of large differential ‘settlements, 


whieh are 


‘TYPICAL TANK BASE DESIGN 


As a result of this experience in the construction of tanks of various sizes” 
on sub-soils of varied characteristics the following general approach has been 
pany’s 8 s standard practice 


evolved. Prior to base construction - is need the c com 
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‘FIG. 11.—SETTLEMENT OF TANK NO. 39 WITH MODIFIED FOUNDATION | 
> 
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TANK FOUNDATIONS 


pro ‘it pon been found to be most economical to investigate the subsoil con-_ 

ditions of all sites to avoid costly failures. With this: data, soil studies are 
_ made to determine the possibility of problems arising from one or more of the 

three critical types of tank foundation settlement resulting from either shear 
or from soilconsolidation:; 


Differential settlements between opposite sides causing the shell to 
form elliptically setting up excessive stresses inthe shell; = | 
2. Differential settlements closely —_—— along the shell causing the shell — 


3. Large differential settlements between the shell and the - 
‘The stability of the entire tank foundation against shear failures 
causing slides, plastic flow and upheaval, or overturning failures is studied. ie 
After determining the soil characteristics and corrective action to be — - 
where necessary a typical base, as shown in Fig. 12, is constructed. gly 4 

With the proper design of connected piping, a typical petroleum tank rel 

signed and constructed in accord with the API Standards can withstand several L 

_ feet of uniform settlement. A large tank canalso operate successfully with the 
center : area settling 2 ft to 3 ft more than the shell because of its flexible bot- _ 
tom and roof. The columns supporting the roof can be adjusted in length from > ; 


time to time to correct excessive deformation in the roof, Experience in the — 


Caripito Farm has amply demonstrated this fact. Nearly every tank 


7 mtx has been opened for inspection after severalyears service has had from 
1 ft. to 2 ft more settlement near the center of the tank than at the shell, > 
+ Site Drainage. —The entire area enclosed within the limits of the firewalls 


built to a minimum of 1 ft above the elevation of the stripped ground surface, : 
This insures drainage of the base even when the area enclosed by the — | 
7 walls is. subjected to heavy rains before the operating personnel can. open 
drain valves. A minimum slope of 1 in. in 10 ft is used to allow for possible 7 
greater settlements of the tank center thanthe tank shell. Both the base thick- 
= ness and finish slope are frequently increased whenit is known that the subsoil : 


conditions may cause considerable settlement G in. or more) as in the — - “4 


: 7 is stripped of grass, roots , and topsoil. The outer edge of the tank base is 


Base Material.—The material specifications for the ‘compacted fill” 
generally limited to Casagrande: classifications GW, GC, GF, SW, SC, and SF or : = 
better, Within the limitations of the material available, uniformity, moisture 
content, and density are closely controlled to insure a minimum of 90% maxi- 
mum density and satisfactory shear strength. The (AASHO) American Assn. of 
State Highway Officials Specification M57-49 for embankments and subgrades 
is followed as a guide for construction. With these simple but flexible specifi- 

cations it has been possible to build successful foundations utilizing local ‘mate-— 
rials in remote areas at minimum cost. __ 
Asphalt Blanket.—A 2 in. mix blanket 
paced over the entire base and finished to within + 1/4 in. of the final grade 7 
This asphalt blanket insures that the finish grade is maintained during con- 
struction of the tank and insulates the tank bottom from the corrosive ef- 
fects of the underlying moist soils during the normal life of the tank. Because 


one o of the most critical ‘Settlement Conditions | of any ta tank is os differential 


to investigate the site by Taking one Or more 
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‘settlement along the tank shell, it is esate that the on pad be well 
x compacted and true to grade along the shell line. The unit pressure directly . 
beneath the tank wall are as high as 130 psi while it is only 22 psi over the id 
interior of the tank bottom in a fully loaded 48 ft high tank. Some tanks’ walls ¥ 
“breathe,” that is, move in or out horizontally as much as 1/2 in. under 
certain loading and temperature conditions. To prevent the outer edge of the - 
4 tank bottom from gradually cutting through the sand asphalt pad, well-graded — 
mixes using RC-2 asphalt are designed. A wide range of mixes may be utilized 
depending upon local materials. However, asphalt content and compaction are > 
both closely | controlled to insure adequate density and stability. — 


ail 1. Careful subsoil exploration, laboratory analysis, and expert design to- 
gether with close control of field construction will provide the most economical — 7 
_ tank foundations, even in areas for whichthereis a great accumulation of data. 
2. Large oil storage tanks can be successfully usedin areas where several 
; feet of settlement is probable, if the settlement is anticipated and provided for. | 
3. The greatest danger to the tank arises from excessive differential 


settlement, especially along the shell. 


_- 4,: The use of concrete ring walls in tank foundations is generally not justi- 
5. The Sai ndahiaaiians of numerous tank foundations in Eastern 
ss has led to the adoption of the compacted granular fill base with 21 in, 
sand-asphalt pad as the most i ania and adaptable for general use. 
ad The writers wish to acknowledge that the information presented in this re- 
port was compiled from the engineering files of the Creole Petroleum Corpora- 
tion Eastern Division.? The original subsoil explorations, laboratory analyses, 
and field supervision involved the efforts of many engineers, however, S. J. 
Mathis, F. ASCE, andC. A, Vickland, M. ASCE, contributed most to the founda- 
tion | engineering examined herein. The work of analyzing the Tank 39 foundation © 
— and the modified foundation was accomplished10- by Ww. Ss. 
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+ “Test Loading of 150,000 Bbl Storage Tanks, Caripito Terminal,” Creole Petrole-_ 

um Corp. Engrs., General Engrg. Report No. 9, January, 1951. 
8 “Report on Tank No. 39, Caripito as ” Creole Petroleum Corp. Engrs. 
General al Engrg. Report No. 15, July, 1952. =< tana 
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‘eee. F, ASCE, and G, A, Oakes, M. ASCE, of the e University of Michigan, 


BASIC DATA 


— _The total test load pressure on original, ground - S. E. side is \ com-— 


‘posed of tank weight of 75 psf, 4 ft earth fill - 420 psf, and 48 ft of water or 
- 995 psf for a total of 3,490 psf, say 3,500 psf. The total test load pressure on : 
_ original ground - N. W. side is composed of the tank weight of 75 psf, 2 ft - 
earth fill or 210 psf, and 48 ft of water or 2, 995 psf, for a total of 3, 280 ‘psf, 

Bearing Capacity Formula. 


in n which qd qd the bearing in pee square foot; No is; a factor 
-_-varying from 5.2 for strip footing to 6.2 for round or square footing - both on 
the ground surface; y denotes density of surcharge, , in pounds per cubic foot; 
d refers: to the depth of surcharge; c is the undrained shear strength; and 
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ing r resistance a circular 


I 1 (Fig. 13). —Assume strip case 


7 


5.24 

or r annular strip 6 ft 


(5.24 x 720) + (-1.3 x 105 f) 

bas x + x pe 

3 - 137 = 3636 pst : 


4 d= (5.24 x 720) + + (2.2 = 400: : 
slightly worse case. 13, b radius of assumed shear 
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_ ‘TANK FOUNDATIONS 


— Case IA.— Check stability to - 6.0ft with base and top stratum 1 replaced wit! with 
material and with 6 ft surcharge 45 ft wide. 


E3490 


factor 


32, 
ad 


= (cylyb + cgLgb) 


b = = 19. ft, = 720 psf, cg = 1860 psf, Ly ‘11. 9 ft, and Lg = 47.8 5 
2 8) = 10,250 psf 
Reduced bearing ng capacity = 0. 86 x 10,250 = 8,840 psf 


of base with 58 ft wide | 


+ + + cy = 1860 pst 


j (1860 x nm. 9+ 1860 x 41. 8) (108 x 6) 


th 

"Reduced bearing capacity = 0.86 x 12, 310 _* 


I A (Fig. 15). .—Check stability from 0 ft to - 73 ft by balancing 
abe ah point Oand reducing bearing value in accord with pannel meth- 


bearing value factor = 


4 2 + colo + » 


— 
= - — 7 
II (Fig. 14).—Check stability from 0 ft to -19 ft by balancing moments 
about point 0 and reduce bearing value in accord with adjusted N, value 
— 
a 
Case II A.—C 


= 


1961 


720 pst c2 = 1860 = 350 psf, 12.8 8 psf, = = 4 and 
qd 73 (720° x ‘12, ms + poe x 20.4 + 350 x 19 196.2) = 3180 psf 


= bearing value reduced = (0.92 x 3180) = * 2926 psf 


Case’ m B.- —Check stability from O ft: to -73 ft with 6 ftx 58 ft wide surcharge 


1860) 0 ft to ft t replaced with material with shear 1860 psf 
1860) 


2 2 + + + cala) + (7 


at 


2 (1860 x 12.8 + 1860 x 20.4 + 350 x 196, 2) + 
6 = 3560 + 398 = 3968 


; a Reduced bearing value = 0.92 x 3968 = = 3650 p 
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NICS AND FOUNDATIONS. DIVISION 


the Ame 


By! Milton E. M. ASCE and Robert. c. 


A type of civil engineering problem involving the flow of groundwater is. 


solved by making use of potential theory, elliptic integrals, the theory of com- -— 
“plex variables, and conformal mapping. _ The results are summarized in chest 
_ form for ready use in analysis, allowing a rapid solution of a problem that re- _ 
quires much time and effort by conventional methods, methods. 


‘Powerful analytical tools been to the solution of a number of 


“Seepage problems (4)(5).3 In the development of these methods, potential 
theory, “elliptic integrals, the theory of complex variables and conformal map-— 
ping have been used to advantage. In spite of their great generality, these 
‘mathematical methods are not widely used in the analysis of structures tor 
which the flow of ground water has a controlling influence on the design. It - | 
the purpose of this paper to show how these methods can be adapted to obtain 
“essentially exact solutions to a variety of civil engineering problems. Although | 

_ an outline of the theoretical ¢ analysis is presented for those who may be inter- —_ 

ested, the primary purpose of this paper is to show how the essential results 
= be summarized in chart form for ready use inanalysis. ~~ 


Note.—Discussion open until March 1, 1962. To extend the closing date one month, a 
: written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Soil Mechanics and Foundations Division, tae 7 
of the American Society of Civil Engineers, Vol. 87, No. SM 5, October, 1961, ~~ 
Assoc. Prof. of Civ. Engrg., Purdue Univ., Lafayette, Ind. 
2 Research Engr. Senior, Materials Research Lab., Kentucky Dept. of Highways, — 
Numerals in n parentheses to © corresponding references in ‘Appendix 
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The problem selected for. is in Fig. ‘The principal 

> factors of interest in design are (a) the quantity of seepage, (b) the pressure 

- in the water at any point, (c) the exit gradient which determines the possibility 
of piping and (d) the average gradient across a critical section that may be 

7 subject to heave . Fig. 1(a) represents atransformed section. The anisotropic 

“a case, in which the | principal axes of permeability | are parallel to the > X- - and y- 
coordinate axes, canbe treated by applying t the appropriate scaling factors (6). aa 


ANALYSIS 
Notation. —The letter symbols adopted for use inthis paper are defined where 
_they first appear, in the illustrations | or in the text, and are arranged alpha- 
r betically, for convenience of reference, in Appendix a 2a 
ee? Essentially, the procedure used in solving such a problem is to transform 
2 the problem from a geometrical domain in which | a solution is sought into one 
for which the solution is known. The mathematical expedient permitting such | 
"7 a transformation is called conformal mapping, which in turn is a } part of the 
of functions of complex variables (2) 
= y) +i w(x, y) = f(z) 
is analytic, then it can be shown that both ¢ and s equation 
and that ¢ = constant and 7 = constant (except at singular points) intersect at 
‘3 right angles in the x- -y plane. It can be shown that in the transformation from 
- “one plane to the other, the relationship w = f(z) may, in general, alter both 


7 the scale and orientation of a Hine, but the ueiae formed by intersecting | lines” 


| between the various on -z-plane and the t- 
7 plane are as shown in Figs. 1(a ) and 1(b). Applying the Schwarz- ‘Christoffel 


Note th that z = 9 when t = 9; thus, N must sisiiiie zero. Eqs. 1 can be expressed = 
as the difference between two integrals in the following manner: 


on 
: The problem then is to find the particular functional relationship that will : Ee 
pee ais . 4 perform the desired transformation. This is generally provided by the Schwarz- aa 
ay Christoffel transformation (2). The Schwarz-Christoffel transformation will Lae 
re _ be used in this problem to map the polygonal regions of the z- or w-planes on Ls 
- 
hose 
— 
— 
— 
2) at » a 


| 


— 
— 
— 
tmz 
(OW PLANE 


Le Letting t t= m7 T, and thus d dt = =m adr, 


and 4 be- 


‘come in 1 Legendre s notation 


“were Fim, 


batting Eqs. 5 and 6 into o Eq. 2b, with 
u-u+ E(u) ]= 


at Point 1, t=, 1 = R(X) * 


At Point 2, z = b+ is, t =a, t/m = o/m > < 1/n, For this case, u 
will be a ‘complex of the form u = =K+iv (1). ‘Thus 


on (iv cn _(iv, m) 


) dn dn (v, 
= dn-1 tee 

in which sn u, enu and dn u are the Jacobian elliptic feiilebiiias (4), and m' is 
the complementary modulus (mt - - m2), It can also be 


E(u) = E(K + iv) = E +i By, 12 


— 
Substituting 7 = sn u= sin =t/m; =am u = sin-! (t/m) 
= F(m, @) - E(m, — 
mf) inteorals of the first and second 
— 
he relationship between 
— 


— 
(sn v', and so forth, designate sn (v, m'), vee) 


+iM 


which K' and are the complete elliptic of the first and second 


kinds with modulus for Point 3 it is found 


=} [ (02. - 1) K+ +iM 1 (02 )=o - id 


_ On this basis, the three constants M, o and m can be determined from 


and o thus can be determined re u 


and substituting from Eq. 14 all. 


— | S 
for Point 2, 
At Point 3, z = b -id, t = 1, t/m = 1/m, u= K+ 
— 
— 13a) a 
| s m is known. 
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FIG. _—DETER MINING THE MODULUS 


— ous 
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— 


(15e) 


- KK' = 1/2), 


substituting the M: and | - 1) into Eq. 7, the follow- 
is obtained the mapping function: 


where the modulus remains 1 to be determined, 
From Eq. 13b, 


which the left part, noted from Eqs. 9 14,. isa of the motdes 


_ and the ratio d/b only. Because of the implicit nature of Eq. 17 a graphical - 
_ procedure was devised whereby the modulus lus could be = for awide range 
of d/b and s/b ratios. Such a plot is given in Fig. 2. rey wr wad aan 


It is noted, in particular, that when s = 0 in Eq. ia. = o, hence Eq. 14 be- - 
“4 


2 k) 
Hence, for s /b = 0, the ‘values of the modulus can be obtained directly from 


is convenient in this problem to assume the as shown in Fig. 1(e 


SM 5 SEEPAGE PROBLEMS 
which, recognizing Legendre’s formula (E'K + EK’ becomes 
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PROBLEMS 
and because w = 0 when ee = 0, N= 0. Making the substitution t=m7T - in Eq. 19, 


the correspondence two planes at Point 1; t = 


we i Q/2 where Q) is the quantity of seepage, then t/m = 1, = 7 


(21) 

wT 

Point 3, t= = 1, w= -k H+ wher 

thm = = “hence 


In addition to the ‘determination of the quantity of seepage the possibility a 


heaving and piping must also be investigated. ‘Heaving of the soil prism 1256, 
Fig. 1(a), 


notation for Eq. 20 


which sin-} (t/m). the a average g gradient along the 
line 62 (Terzaghi’s method) (7), the gradient along the line 12 will be obtained. 
The factor of safety with respect to heaving determined on this basis will a 
ways be slightly on the safe side. 

At Point 2, w= +i Q/2, t = o and t/m = 1. Because F( 


and for the ne gradient along th piling, 


We 


— a 
— . 
— 
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_SEEPAGE PROBLEMS = 101 
in which v = dn-1 i sidilin m' ), Eq. 9, with o as determined from Eq. 14, Graph- 
ical solutions for { ip s/H of Eq. 26 are given on in Fig. 5 for various rettes of d/o. 
modulus can be obtainedfrom Fig. 2, 
Considering the exit gradient (Points 1 and 5), 


126 at az 


i 


3 o2 - 1 E 
plot of Eq. 28 for ie 8/H is given in Fig. 
the the information necessary | to effect a is now in graph. 


| 


illustrate the use « charts in solving the con-— 
Rey illustrated in Fig. 7 will be analysed. For example, when s = 5 ft, b/s 
= 1.00 and d/b = 2, from Fig. 2, m = 0.108. “Entering. Fig. 4 with this value for 
: m, Q/kH = 0.90, from > which the ssinia of seepage is found tobe 13.5 cu ft per 


‘From Fig. 5 (with m = 0.108, ip s/H = 0.54, that is, ip = 1.62, The factor 


of safety against heaving nemealen a —— along the inside of the sheet piles) 


_ From Fig. 6, ie 8/H is found to be 0.40, whence ie = 1.20, and ‘eine 
safety (using the exit gradient) is equalto 0.72. 
This procedure can be repeated for various depths of penetration of the 
_ sheet piling, and plots of the quantity of seepage and the factor of safety with | 7” 
respect to heaving versus the depth of penetration can be obtained as shown in 
a me The methods used in solving the illustrative problem (Figs. 1 and 7) can be 
extended to other situations, such as the one illustrated in Fig. It is 


the ‘same mapping from the z-plane to the t- -plane is as given in ay 1. How. - 
ever, it is noted that the w-plane is different. Fig. 9(c) indicates. that the 

quantity of is infinite 1(c) indicates finite quantity of 
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FA CTORS OF SAFETY AND QUANTITY OF 
SEEPAGE AS A FUNCTION OF DEPTH OF 


SHEET PILE PENETRATION 


“25.1. = 


os 


: — 


W-PLANE 


RESSED WEIR WITH SHEET PILES 

OF Verve 


> 


spa cos-1 


and hence the in the sted) 
and 
For the exit gradient (at Point 3), 


m2 _ 


Noting that t = at exit and for or (02 - - 1 from Eq. 15b, 
the exit gradient is | given minh 


where the modulus can be obtained from Fig. 2. _ 


_ Consider now the gradient along line 23. At Point hea , from Eq. 29b, 


T 
is and for the gradient along the piling, 


Se 


ere m is ‘obtained from Fig. 2: and o from Eq. 


ae It is — that this paper has illustrated how the powerful analytical tools © - 
available to civil engineers canbe applied to obtain “design” curves for various | 
problems. These methods have been used in this paper to develop a group of 
‘4 design curves that can be quickly and easily used to aid in the solution of a 


type of civil: problem involving the flow of groundwater. 


owe following naman have been adopted for use in n this paper: 


= 


— 

‘ 
— 
— 
— 


me weir is or depth of excavation between sheet pil- 


elliptic of the | second kind with modulus m; 


= complete elliptic integral of ‘the second kind with modulus im; 


es = elliptic integral of the ‘second kind; ae | 


factor o of safety; 


mathematical operator, the imaginary unit; 
= critical radient; vr, 
= exit gradient; 


complete elliptic integral ‘of the first kind with modulus 
‘* complete elliptic integral of the first kind with modulus ms 


Ver 
= co- modulus; ieee 


= constant of integration; _ 


= pressure in the 


- 
depth of sheet pile | penetration; 
variable, a point in the t-plane; 


‘complex a in the w w-plane; 


complex: variable, a point in the z- plane; 


unit t weight of water; ‘ 


— 
— > 
— i — 
— 

= gradient along the sheet piling; om 

 «K 

— 

> 


SEEPAGE PROBLEMS» 
= coordinate of a 2 point ‘the t- 
= dummy variable; 


abscissa of point w; and 
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eli DRIVING semen BY THE WAVE EQUAT 


by 


aes of computation that would not materially affect the: final result but are 


It is ‘characteristic of approximate solutions that there is iain a choice ° 

4 as to how the approximation is to be made. Soderberg has developed a method 

- somewhat different from that of the writer, and has set sort Se ens essentials in 

his analysis, He should have full credit for this development. 

—— Soderberg computes only the displacements. Forces or other information 

_ that may be wanted must be obtained i by means of a supplementary computation, — 
therefore in many cases the time | saving might be much less than Soderberg > 
has indicated. In any event, this is a subject for study primarily by — | 


_ticians and computer experts as it does not affect the civil engineering aspects 
a Gates examines the question « of how the bearing c capacity of piles should “7 
determined. For this purpose Gates opposes the use of the wave equation, but . 

- advocates a pile driving formula (Eq. 37) which he has devised. The writer 

_ presented this general subject briefly in the paper under the headings “Intro- | 

_ duction” and “Recommendations and Comments .” The words were carefully 
chosen, and there is no wish to addto them herein or to amend them in any way. — 

4 In Gates discussion there are three statements that call for answers -" 

 e The discussion states that previous investigators “unanimously ca tioned 
against the use of the wave equation asa means of determining bearing capac- 
” No specific references are givento support this statement, and the writer 
has been unsuccessful in his attempt to confirm it. He did, however, find ~ 

_ statement exactly to the contrary in the writings of A. B _ Cummings43 whom > 

both Gates and the writer have referred to as an authority. This same - state- : 
am also appears in the last part of a paragraph on page 13 of Gates’ refer- 

ence No. 25. Cummings first refers to the wave equation developments of St. 
mt Boussinesq, Isaacs, and the British Building Research Board (Gates’ 

references Numbers 11, 12, 17, and 18) andthen writes as follows: “However, . 
there is a considerable amount of field evidence available which shows that 
the stress transmission characteristics of a pile are of great importance not 


ia August, 1960, by E. A. L. Smith (Proc. Paper 2574). ee ee 


42 Chf. Mech. Engr., retired, Raymond Internatl. Inc., New York, N.Y. 
_ 43 Discussion by A. E. Cummings of “Pile Driving Formulas, ” by Progress Report — 
of the Committee on the sens Value of Pile Foundations, Proceedings, ane, ™ 68, 
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only in determining its behavior during driving but also with : respect to its sub- 
sequent ability to carry static load.” Cummings’ paragraph c concludes: 
and promising field for investigation.” 
_ B. It is suggested in the discussion that the writer “should make a definitive h 
concerning the limits of application ofhis presentation.” writer = 
has no intention of doing this because he believes that there is no good reason 
for limitations beyond those already stated inthe paper. 
ae _C, It is stated in the discussion that the writer is guilty of “qualification of 
be his previous references to bearing capacity by mentioning in this paper soils 
_ which may either relax or set-up.” However, the writer had previously men- | 
wae 7 tioned these particular soil characteristics in the second paragraph of Gates’ a 
a reference No. 16 published in 1957. No such 1 “qualification” has occurred. 


7 Gates has implied that the wave equation and a dynamic pile formula are 


* : tion my ’ therefore be cited against it. The wave equation is a mathematical © 
"method. ? Dynamic pile formulas appear to be ‘mathematical, but their mathe- 
matics is defective as has been carefully explained by Cummings3, 43 intwo — 

: ow Passing completely from the subject of bearing capacity, the analysis next = 
7 ‘ analysis makes statements which, if taken together, leave the reader with the : 

_ impression that the writer’s method is quite inaccurate. These statements are ; 


1: The discussion lists six necessary in order to apply 
the wave equation to pile driving. These have been taken from a list!9 that 

a applied specifically to the wave equation pioneering done almost 100 yr ago 


by st Venant!1 and 12 This list is now (1961) out-of-date. Caly 


caniiiaen but it is not necessary because side friction has been oie 
‘provided in the writer’s method. , Assumption No, 2 involves negligible error 


if if the ram is: short and stubby, but a comparatively long r ram may be see 


= 
‘been provided for. No. 6 is completely unnecessary. 
_ _ No, 2: Immediately following the list of assumptions it is stated that. thatthe 
i writer’ s method takes no account of propagation losses in the pile. The writer 
y has presented ‘Eq. 27 for this purpose. However, in piles of ordinary length the 
‘Ne. 3: Near the end of the discussion, ‘the. compressive forces computed in 
pi the writer’s Illustrative Example are compared unfavorably wi with those obtained a 
A id by strain gage measurement in certain field tests. 30 The discussion emphasizes i. 
ea ‘that the computed forces increased rapidly towards the tip of the pile, and of 
states “This conclusion is diametrically opposed to the results obtained from 
field tests.” Gates has | ignored the fact that in the Illustrative Example, side — 
friction was completely absent, whereas the field measurements were made on 
piles in which side friction predominated. There is no reason whatever to 
expect the forces to be distributed in the same manner in these two radically 
zy different cases. The discussion also emphasizes that the computed forces are 
_ larger than the measured forces. Here again a vital fact is ignored, namely 
a that in the two cas cases the pile weight per foot is different. The heavier and a Tia 
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DISCUSSION 


stiffer pile used in the Example will inevitably higher im- 
pact forces if the same hammer is used. A comparison of this kind should not 
have been made because the only major way in which the two cases resemble 
each other is that the same hammer was used. 
~~ _ No. 4: The discussion states that inthe writer’s method there is a “need to 

- 4 preassign not only the factor R,, but also the penetration at which this value 

_ occurs, ” As a matter of fact the penetration “s” is obtained by making the ; 

"wave equation computation. It is never under any circumstances | “preassigned. oil 
This appears in the paper as part of the first paragraph under the heading = 
“Outline of Numerical Method” and also as part of “Conclusion No, 1.” It also 

appears in “Step No. 9” of the Illustrative Example it the 


namely: should» Ry b e 200,000 Ib as given in the paper, or should Ry be a 

| 000 Ib or 405,000 Ib as suggested inthe discussion. = = 
The discussion states: “It can be shown that in the Illustrative Example . 

=z by the author, at me end of set 


Jonas has written a thoughtful discussion which raises s an question, 


as can ‘readily be seen referring to Fig. 12. At this instant 
Fu (the: internal force at the point of the pile) does not have its maximum = 
value" of 405,000, as quite | reasonably assumed by — but has value of 


statenee. ) This fact can be determined by reading mse from the curves of 
_ Fig. 12. Thus in time interval 55 Dj, = 0.202 ‘in., , Dyg = 0.302 in. and Dj; = a 
0. 353 in. Subtracting Dp from DjQ gives 0.1 in. as ‘the measure of the elastic - 
compression. Multiplying 0.1 in. by the spring constant of the ground» 
Kp (= Kj2 = 2,000 000) gives the value of Ry as 200,000 lb. Subtracting Diz 


point of the pile’ 3, ,900,000) gives Fi, as 200,000 Ib also. Thus equili- 
s The maximum value of F;,; (405,000 Ib) occurs not “at the end of permanent 5 
_ set” but in time interval 34 when the point of the pile is forcing its way through | 
ot the ground with the rather high velocity of 6.4 fps. An analysis similar to that 
. given in the preceding paragraph but involving resort to the printed computer 
output, shows that in time interval 34 the values 
Ultimate ground resistance Ry... .. 200, 000 Ib 
‘This total is equal to ‘the maximum internal force in th Fit. 
es Therefore, it may be concluded that the correct value of R,, is 200,000 Ib 
that no modification of Eq. iS and 16 is required. 
a) Eremin’ s discussion raises a question regarding the velocity of sc sound 
im fl (or stress) in a pile. None of the - equations used in the writer’ s method in- 
4 cludes, as a term, the velocity of sound. Nevertheless, as the computation is ; 
_ performed, the stress wave progresses down the pile with approximately the 
we of sound, as it should, Furthermore, if side resistance is included, _ the 
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stress wave will be slowed down alata Eremin states that it one 
be. This is an important proof of the correctness of the method. a 
‘The values of 0.05 and 0. 15 proposed for the damping constant J' and J are 1 
the best now (1961) available. When experimentation or other means produces» a 
suggessts that “tip bearing stresses,” “buckling stresses,” and “elas- 
= tie properties of soil” have an effect onthe damping forces. The effects of the 7 a 
ia first and the last of these are considered in the paper and result in Eqs. 15 and - 
16, and Eqs. 28 and 29. After the pile enters the ground there is little — 4 


— to buckle if the pile is reasonably straight, therefore the effect of buckling an 4 
stresses has not been considered in the writer’s paper. we 


_ Eremin refers to data obtained from driving piles by vibration, No mention 
‘was made of such data in the paper. 
 Cornfield’s discussion is more than welcome. The English were pioneers in | 

_ applying the wave equation to pile driving.5 However, the fact that —— = 
- computers were not available in 1938 when tl their major paper? was written, - 


. The first paragraph of Cornfield’, s discussion contains the essence of a lot 
of practical pile driving eupantonce, ont is worthy of careful 
parameter investigation is suggested in connection with the writer’ 
method. Samson of Texas A & M College is at present (1961) planning such i 
_ Cornfield draws certain tentative conclusions, partly from wave equation © : : 
_ data previously published35,36 and partly from test load results. ‘Conclusions 
such as these are valuable if applied within specific limits but may become — 
inaccurate beyond t these limits. The wave equation can be helpful in Gotermining i 
. The curves in references 35 and 36 were prepared before viscous eile 
_ was introduced into the wave equation computation, Inclusion of viscous damp- — 
“S Chellis and Zaskey have presented an extensive and important discussion. - 
The writer commends it highly. His only comments are by way of adding — 
As a means of saving space, tabulation has been found preferable to curve > 
- drawing for recording the results of a considerable number of wave equation © 
| The wave equation curves of Figs. 15 and 16 are based on the assumption © 
=J' = 0, whereas it is recommended in the writer’s paper that until 
more accurate values are available, J should be given the value of 0.15 and > . 7 
#£ the value of 0.05. Consequently these curves (as noted in the discussion) _ 
oy yield values that are less conservative than those that would be obtained by | 
7 The programming of a computer at Texas A & M College, College Station, 
7 ‘Texas, for wave equation computations, as mentioned in the discussion, is. 
under Charles H. Samson, Jr, A. M. ASCE. Use of the computer and program, 
which contains certain refinements introduced by Samson, willbe made avail- 
a able to anyone. The project is being sponsored by the Texas Highway Depart- 


ment and is under Randle B. sae, #. ASCE, Bridge Engineer, _— 
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"DISCUSSION 


In — writer w wishes to thank all who presented discussions, 


< ‘They have thrown light on many aspects of the that t were e not 


38, under DIA DIAGRAMMATIC REPRESENTATION the v word 
7 “ean” should read “cap”; page 44, line 14, “Eq. 24” should read “Eq. 27”; 
page 47, line 3, “Q one” should ‘read “Q”; page 50, three lines from bottom, 
“Slide” should read “Side”; page 61, line 5, “Eq. 24” should read “Eq. 27”; 
 ~page 61, line 7, “25, and 26” should read “28 and 29”; page 61, four lines from 
bottom, “Eq. 24” should read “Eq. 27”; andpage 61, next to last line, “Eq. 24” 
Discussion by Chellis and Zaskey. —Page 67, line 12, “impluse” should 
“impulse”; page 67 , line” 31, should read “or”; page 67, line 36, 
“lone” should read “long”; page 69, line 16, “effects” should read “affects”; ; 
age 69, — 39, D Delete extra letter 6" ; and d page 74, line 2 18, “easy” should 7 


“each.” 
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H, Grice 


H. GRICE.29—A Penetration Ratio 
of grouting fine fissures is presented. The writer is particularly interested in 


m how this ratio was obtained ashe considers that it is a fundamental character- 

Ny 4 istic of grout mixes. Its measurement shouldbe an essential part of a standard- 
_ ized procedure for use in the design and quality control of grout mixes. — 
a “; penetration test would be of particular value in the study of lubricated — 

grouting and its relationship with dispersant treated grouts. Work on these 
4 lines would be likely to modify the groutability criteria proposed in this paper. oi 


Mee Test standardization is particularly important as mixing procedures, injec- 
2 ‘tion pressures and rates are critical in the design and handling of the highly 
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RICHARD R. SCHINDLER, 12 M. ASCE. elim-_ 
of foundation vibrations by injecting chemical grout to alter the elastic 
properties of the soil and the natural frequencies of vibration of the foundation-_ 
; _ soil system, Such an approach may be of interest also to the structural engineer 
- concerned with the prevention of damage due to seismic disturbances or nuclear 
_ weapons e energy releases. This problem differs in being one of transient = 
= in the magnitudes : of energies and weights involved. But the conclusions of 7 : 
Gnaedinger and the data of the Symposium appear pertinent. Also, the type of 
soils correctly associated with earthquake damage a are, fortuitously and fortu- 
_ This discussion attempts to evaluate qualitatively the effects of grouting dl 


4 of consolidation grouting, could contribute to forestalling 


earthquake damage or to reducing the need for ‘structural provisions ag —- 
Intensity of Ground Motion in Relation to Soil Properties. — The change 
energy of seismic waves passing from one medium to another is still a matter 
of considerable study. But B. Gutenberg!3 cites extensive data indicating that 
the intensities perceived by instruments situated on rock are less than those . 
recorded by instruments on alluvium. (This applies primarily to later waves, — ‘4 
; the initial _ shock being usually recorded with relatively little difference.) He — 
that, “Ground effects: may produce appreciable differences in 


to double in intensity, as shown theoretically and substantiated with field meas- ; 
urements by G. W. Housner, 14M. ASCE. It is not inconceivable that passing a 
Ss alluvium to rock the intensity would decrease. But the degree of similarity — _ 

_ between a grouted soil and bedrock needs tobe determined. = 


The effect of grouting « on the intensity of seismic oscillations should be con- sof 
sidered negligible until soundly established. Yet it is likely to be beneficial ed 


‘rather t than de 


— @ April 1961, by J. P, Gnaedinger (Proc, Paper 2793). 


Paul J. Tolen, Structural Engr., Los Angeles, Calif, 


_-:13 “Effects of Ground on Earthquake Motion,” by B. Gutenberg, Bulletin, Seismolog- — 


ical Soc. of America, Vol, 47, No, 3, July, 1957, 
14 “Geotechnical Problems of Destructive Earthquakes, " by G. W. Housner, Geo- 
technique, Vol. 4, No.4, December, 1954, 
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Soil Interaction. —A building on agrouted foundation, with in- 
creased elasticity will vibrate with a shorter frequency than on the original 
5 soil. This affects the stress conditions inthe structure adversely. The gal 
: _—-panying smaller rocking displacements tend to aid the stress conditions, but, — 
% Bn simultaneously, may reduce desirable damping effects. Energy dissipation by 
_ stress waves produced in the ground by oscillating buildings has not been ovalu- 
_ Analyzing the effect of foundation conditions over a wide range in their inter- : 
. action with five, ten and fifteen story buildings, Merritt and Housner15 - 7 
_ Clude that rocking motion is not a significant factor for most ordinary struc- é 7 
ae, tures, frequencies being altered by less than 10% over a range of all foundation — = 
- conditions to be anticipated in practice. The extraordinary structures are only 
those which are both narrow andon soft soilh 
‘The effect of grouting on the ol -aalieebtrwteitis interaction, which is so 


important for elimination of vibration, appears, consequently, not to be signif- — 


damage is not only the intensity of the shock per se , but also the differential — 


f et on rock or hard ground. Such studies also indicate that the cause of a : 


_ Although many earthquake damage reports come to this conclusion (Mexico 

1957, Japan 1948, Tehacha 1 Thet and others), aparticularily striking example 

‘given by V. A. Murphy:! 

= ly on an85 ft deep alluvium and partly o on a 100. ft high limestone hill. ‘During | a is 

a 1931 earthquake the buildings located on the hill experienced no serious founda- a a 
failures” only light damage to the structures. ‘The buildings on the 


and damage to the superstructures. The lesser damage on the hill could well 
be due to the lesser shock intensity transmitted by rock. However, emphasis 
is given the observation that the failures on the alluvium are, judging by their 
7 nature, attributable to foundation settlement rather than to building motion. 
‘gaat That densification by grouting, in effect analogous to compaction, decreases — 
4 4 the settlement potential of a soil is almost self-evident, and is illustrated by 
Bruce E. Clark, 17 EP 
Where unfavorable soil conditions obtain, reduction in differential settlement 
* resulting from consolidation induced by seismic vibrations and building motion © 
may be achieved by grouting the soil. The engineer has to balance the factors 
a pertaining to the grouting operation with those connected with making struc- 
tural provisions against the differential settlement. 
: ‘The three aspects presented indicate that in aseismic design grouting may | 
- find another field of application. But the solution of soil mechanics a : 


ce connected with either steady state o or r transient motion are still a a 


«15 “Effect of Foundation Compliance on Earthquake Stresses in Multistory Buildings,” — 
Merritt and Soc, of America, Vol, 44, No. 4, Octo- 


Consolidation —Studies of earthquake damage have demonstrated tha 
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a 
“The Effect of Ground Characteristics on the Aseismic Design of Structures,” by 
et 7 V. A. Murphy, Proceedings, 2nd World Conf, on Earthquake Engrg., Japan, a ia 
ii __ 17 “Symposium on Grouting: Grouting at Ft, Campbell Theatre,” by Bruce E. Clark, i ae 
_ Proceedings, ASCE, Vol. 87, No.SM2,April,1961,. 


DISCUSSION 
y many investigators. 18,19,20 a systematic approach based on Quinlan’s 
_ elastic Theory of Soil Mechanics?! appears quite promising. It could lead to. 7 


vibrating foundation-soil systems. By an extension this approach could perhaps 7 
also be made applicable to the determination of natural frequencies of impul- $ 

4 sively disturbed foundation-superstructure systems, sought by engineers con- 

ub. fronted with earthequake or nuclear weapons hazards. Even without a rigorous - 
solution, the experience with and research on grouting published in this Sympo- _ 


demonstrates its value in vibration of foundations 


“Bibliography of Effects of Soil Conditions on Earthquake Damage, Duke, 
Earthquake Engrg. Research Inst.,1958. j= 
“Symposium on Dynamic Testing of Soils,” Amer, Soc, for Testing Materials, 
20 “Proceedings of the Second World Conference on Earthquake Engineering,” Science _ 
f “The Elastic Theory of Soil Dynamics,” ‘by P. M. Quinlan, Symposium « on Dynamic 
‘Testing of Soils, Amer. Soc. for Testing Materials, 1953, 4 
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UNCERTAINTIES IN INE EARTH DA DAM DESIGN4 


who had the opportunity to examine | the stability of earth dams built in —s 

_ several centuries ago. 8 The added interest is from the fact that there have been | 
some catastrophic failures of some recently built earthdams and one example = 
is briefly discussed in the following. 
y "The recently completed Panshet earth dam which was expected to augment : 
a: drinking water supply of the Poona City, has died a pana death, 760 ft = 


rainfall in the catchment area of the Mutha river is somewhere around 200 in, ; 4 
and the earlier expressed doubts (expressed in 1957 before starting of con- | = 


_ struction) whether high earth dams in this region which is particularly defi- 
~ cient in suitable earth for the core as well as the casing, would be safe in > 
‘heavy rainfall , came true, The floods caused by the breach in the Panshet dam > 
meant the natural death, 100 ft top to bottom breach, of the Khadakvasala 
- dam which is 12 miles downstream and which had served well for — 
a several decades. The breaches" in the dams were responsible for the sudden = 


upsurge of water in the Mutha > rendering about 60,000 people homeless with © 
‘loss of life placed around 200. This only explains the one other uncertainity— 
danger from overtopping during monsoon—which is inherent in earth-dam 
design. One of the logical developments in earth-dam construction of special 
significance to such regions will | ‘be towards self-spillway rockfill type.9 
_ _The final sections adopted for the embankment of Terminus Dam and New 
» - Hogan are no doubt among the best possible solutions. The writer, however, — 
would like to refer to two factors, viz. the conservative nature of estimating — 
Bendel properties of soil and the very low values of safety factors allowed _ 
g the design, which perhaps have a balancing affect on one another to have 
is reported that cohesion has been neglected on the basis that it is not bea 


4 


considered to be a unique property of the soil. . . .and that the minimum value — 

of ¢ from the unconsolidated-undrained tests has been arbitrarily selected as 
e 11°, It is tobe pointed out that this muchavowed conservation is unrealistic -_— 
mt contrary to reality. The C'= Ohypothesis in long term satability analysis10 


& June 1961, by John M. Bird (Proc, Paper 2832), _ wera 
7 a Dept. of Building Science, Univ. of Liverpool, United iKingdom, 
oa 8 “Stability of Earth Dams Built in the Past in India in the Light of Modern Develop- 
ments in Soil Mechanics,” by B, V, Ranganatham, Transactions, C.B.I.P., India Research 
_ 9 “Flow of Water Through Rockfill and its Application t to the Design of Dams,” ” by . 
x. Wilkins, Proceedings, 2nd Australian-New Zealand Conference, S, M. and F.E., 1956. 
_ 10 “Effective Stress Method of Stability Analysis to Determine the Cohesion Intercept 
- re at Long Term Failures,” by D, J, Henkel, Proceedings, 4th Internatl. Com., 
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is an emperical one derived from observations of _— failures in overcon- 
solidated fissured clays with load release (as in cuttings) and with access of 
adequate free water. The value of ¢" ‘used with this hypothesis is that deter- 
.. mined from consolidated-drained tests. In addition, anearthdam, besides being 
far from the conditions warranted for the validity of the preceding hypothesis, — 
4 oo with time owing to (1) dissipation of any excess construction 
pore pressures (thereby there is a gain in strength, if any) and (2) the in- 
ee chemical and structural bond perhaps with the deposition of fine 
It is also reported that the design value of the angle of internal friction for 
: 2 the rockfill has been reduced to 40° asameasure of safety against the develop- 
ment of potential sliding planes in the rock shells which could result through 
_ placement of thin continuous layers of the mica-like particles. This approach — 
_ seems to be purely a speculative concept unrelated to any experimental evi- 
3? dence or field observation. The writer would like to suggest the investigation © 5 
Of shearing strength of saturated meta-sandstone with flat mica-like grains by 
-: "conducting field shear tests!1 with the aidof square steel frames with an area 
of about 1.0 sq m or more and a height of about0.3m. || 
_ With regard to the consideration of field methods to minimize nonuniform- 
5 ities, a minor observation will be made of the successful blending of materials a 
in a recently constructed earth dam, by simultaneous supply from more than © 
Si - one borrow area at a predetermined proportion, The borrow areas were demar- 
cated dependent on the uniformity (statistical) of grading of the materials. Int ; 
yi addition spot checks were kept continually during the progress of construction. a 
ca, Thus, in the writer’s opinion, the extensive survey of grading of borrow mate-— 
_ vials now reported in the paper may be usefully followed up with | obtaining the — 
- - actual grading of the compacted fill in the impervious core of the New Hogan : 
Dam at least for reach. Such data serve useful for 


11 “Field Tests on Soils and Rocks, by A. , Dvorak and Peter, ‘Proceedings, sth 
Internati, Conf., 8. M. . and F. E., Vol. 
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